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We re-examine the Lemaitre-Tolman-Bondi (LTB) solutions with a dust source admitting symme- 
try centers. We consider as free parameters of the solutions the initial value functions: Yi, pi, ^IZi, 
obtained by restricting the curvature radius, Y = ^/gee, the rest mass density, p, and the 3- 
dimensional Ricci scalar of the rest frames, ^1Z, to an arbitrary regular Cauchy hypersurface, 
%, marked by constant cosmic time (t — ti). Using Y% to fix the radial coordinate and the topology 
(homeomorphic class) of Ti, and scaling the time evolution in terms of an adimensional scale factor 
y = Y/Yi, we show that the dynamics, regularity conditions and geometric features of the models 
are determined by pi, ^TZi and by suitably constructed volume averages and contrast functions 
expressible in terms of invariant scalars defined in %. These quantities lead to a straightforward 
characterization of initial conditions in terms of the nature of the inhomogeneity of %, as density 
and/or curvature overdensities ("lumps") and underdensities ("voids") around a symmetry center. 
In general, only models with initial density and curvature lumps evolve without shell crossing singu- 
larities, though special classes of initial conditions, associated with a simmultaneous big bang, allow 
for a regular evolution for initial density and curvature voids. Specific restrictions are found so that 
a regular evolution for t > ti is possible for initial voids. A step-by-step guideline is provided for 
using the new variables in the construction of LTB models and for plotting all relevant quantities. 



I. INTRODUCTION 



Lemaitre-Tolman-Bondi metrics with a dust source and a comoving and geodesic 4- velocity constitute a well known 
class of exact solutions of Einstein's field equations |21| ^2j. Under the prevailing conditions of a present day, 



standard, matter domminated universe, a dust momentum-energy tensor is a good approximation to non-relativistic 
matter sources. This fact, together with their mathematical simplicity, makes LTB solutions adequate models of 
inhomogeneous configurations of non-relativistic and collisionless matter (either barions or CDM). The relevance of 
such models is evident as they complement the usual perturbative approach by allowing one to study a case of exact 
non-linear evolution of inhomogeneities [[l0| [^| [^2| [^9). Also, recent literature M Gj| illustrates how a nearly 
isotropic CMB might not rule out an inhomogeneous universe compatible with current CBR observations. 

The process of integration of the field equations for an LTB metric with a dust source in a comoving frame is very 
straightforward. This process, reviewed in section II, leads to a Friedmann-like dynamical equation for the curvatue 
radius Y = ^/g$e, containing two arbitrary functions of the radial coordinate: M(r) and E{r). The integration of the 
dynamical equation yields a third arbitrary function, t bb (r) . The first two of these functions (M and E) are usualy 
identified as the relativistic analogues of the newtonian mass and the local energy associated with the comoving dust 
shells, the third function (t bb ) can be interpreted as the "bang time", marking the proper time for the initial (or 
final) singularity for each comoving oberver. The functions M, E, t bb are then the free parameters of the solutions, 
though any one of them can always be used to fix the radial coordinate, thus there are realy only two free functions. 
Regularity conditions, related to well behaved symmetry centers, abscence of shell crossing singularities and surface 
layers, can be given as restrictions on these functions and their radial gradients. These functions are usualy selected 
by means of convenient ansatzes based on the interpretations described above and/or the compatibility with regularity 
conditions or observational criteria. Once the free parameters have been selected we have a completely determined 
LTB model. 

Since the standard free functions lead to a consistent description of LTB models, these free functions are used in most 
literature dealing with these solutions (see p^] for a comprehensive and autoritative review of this literature). There 
are papers considering specific modifications of the standard free functions in order to formulate initial conditions 
for studying the censorship of singularities poj , or aiming at an "intial value" formulation of LTB models ||]. We 
believe that new variables, defined along an arbitrary and regular "initial" hypersurface t = ti, lead to a more intuitive 
understanding of LTB solutions. Therefore, we propose in this paper using a new set of basic free functions along the 
following steps 

• Consider the rest mass density, p, the 3-dimensional Ricci scalar, ^7Z, and the curvature radius, Y, all of them 
evaluated along an arbitrary Cauchy hypersurface, %, marked by constant cosmic time t = ti. This leads to: 
Pi(r), ( 3 >lZi(r) and Yi(r), where the subindex i indicates evaluation along at t = ti (we shall use this convention 
henceforth) . 
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• Rescale Y with Yi, leading to an adimensional scale factor y = Y/Y. L , so that we can distinguish between the 
regular vanishing of Y (at a symmetry center, now given as Yi = 0) and a central singularity (now characterized 
as y = 0). Likewise, we can distinguish between the two different situations when Y' — 0: regular vanishing 
or surface layers if Y( — for a single value of r (thus restricting gradients of pi and ^IZi) and shell crossings 
when T = (Y'/Y)/(Y[/Yi) = 0. 

• Among the three free functions, Pi(r), ^lZi(r) and Yi(r), we fix the radial coordinate (as well as the topology 
of %) by a convenient selection of Yi, using the remaining two functions as initial value functions. 

• From the latter functions we can construct suitable volume averages and contrast functions, all of which ex- 
pressible in terms of invariant scalars defined in %. These auxiliary functions provide a clear cut and precise 
characterization of the nature of the initial inhomogeneity around a symmetry center, as overdensities ("lumps") 
or underdensities ("voids") of pi and ^TZi. 

The old functions M, E, t bb and their gradients can always be recovered and expressed in terms of the above mentioned 
volume averages and contrast functions. Hence, all known regularity conditions given in the old variables can easily 
be recast in terms of the new ones. In particular, the known conditions for avoiding shell crossing singularities |lq ] 
(see also [[t6) and [0) have a more appealing and elegant form with the new variables, since with the old variables the 
signs of the gradients M', E' , t' do not distinguish manifestly between lumps or voids. In terms of the new variables, 
it is evident that, in general, shell crossings are avoided only for lumps of pi and ^IZi- However, some particular 
cases of initial conditions, such as a "simmultaneous bang time" (t' = 0), do allow for some initial conditions given 
as voids to evolve without shell crossings. If we consider an evolution in the time range t > ti, then (under specific 
restrictions) initial voids can also evolve without shell crossings. All these features could also emerge with the old free 
functions, but the new variables provide a more straightforward and intuitive characterization of the effect of initial 
conditions in the regular evolution of LTB models. 

The new approach to LTB dust solutions presented in this paper accounts, not only for an initial value treatment 
on a regular Cauchy hypersurface %, but to a new formalism based on rescaling the evolution of the models to 
variables defined in this hypersurface. Some aspects of this formalism (average and contrast functions) have been 
developed elsewhere (see and j28| for applications of LTB metrics with viscous fluid sources) . Variables similar to 
the ones introduced here have been defined in [|j , [^0| , ||] , and [ fl [ (the latter in connection with Szekeres dust 
solutions, see equations 2.4.8 to 2.4.14 of ]2^]). Mena and Tavakol |llf| have claimed that some of these variables are 
not covariant, thus suggesting covariant expressions for density contrasts (see also Q). Although we show that the 
new variables can be expressed in terms of invariant scalars, and so are coordinate independent and can be computed 
for any coordinate system, we do not claim that the density and curvature averages and contrast functions introduced 
here have a clear "covariant interpretation" , or that they are unique. However, what we can certainly claim is that 
these new variables emerge naturaly from the field equations, provide an adequate initial value treatment and yield a 
clear and concise description of regularity conditions. 

The paper is organized as follows: section [n] reviews the basic expressions (metric, evolution equation and its 



solutions) characteristic of LTB models in the usual variables. The new variables are introduced in section III 
leading to new forms for the metric, the rest mass density and the solutions of the evolution equation. Volume 
averages and contrast functions along % are defined in section IV . Section [y] examines in detail generic prop ertie s 



of the new variables: regularity, symmetry centers and geometric features of the hypersurface % (subsection VA), 



characteriz ation of singularities (subsection V B), choice of radial coordinate and hom como rphic class (topology) of Ti 



(subsection V C ), expr essing the new variables in terms of invariant scalars (subsection VP ), FLRW and Schwarzschild 
(subsection limits VE). In section [Vl| we look at the characterization of the inhomogneity of % in ter ms d ensity and 



3-curvature "lumps" and "voids" and its relation with the contrast and average functions. In section VII we use the 
new variables in order to look at the conditions to avoid shell crossings, big bang times, regularity of the selected 
hypersurfac e % a nd simmultaneous big bang times. A summary of the no-shell-crossing conditions is given in Table 



1. Section VIII provides simple functional ansatzes for the new variables and their auxiliary quantities (volume 
averages, contrast functions, as well as the old variables and the big bang times). A simple step-by-step guideline 
is provided for the construction of LTB models and for plotting all relevant quantities derived from these ansatzes. 
The characterization of intial conditions as "lumps" and "voids" , as well as the fulfilment of the no-shell-crossing 
conditions is illustrated and tested graphicaly in 12 figures that complement the information provided in the text. 
Finaly, the Appendix provides a useful summary (in terms of the usual variables) of the generic geometric features of 
the solutions (regularity conditions, singularities, symmetry centers, surface layers). 

Some important topics have been left out in this paper: time evolution of lumps and voids by generalizing volume 
averages and contrast functions to arbitrary times (t ^ ti), models without centers, models with a mixed dynamics 
(elliptic and hyperbolic, etc), matchings of various LTB models, etc. The examination of these topics with the 
formalism presented here is currently in progress p9[ . 
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II. LTB SOLUTIONS, THE USUAL APPROACH. 



Lemaitre-Tolman-Bondi (LTB) solutions are characterized by the line element 

y/2 

ds 2 = -c 2 dt 2 + - - dr 2 + Y 2 [d6 2 +sm 2 (0)d4> 2 ] (1) 

where Y = Y(t,r), E = E(r) and a prime denotes derivative with respect to r. The usual dust source associated 
with (0) is 

Tab a b /r>\ 

= p U U , (2) 

where the 4- velocity is comoving u a = cSf. and Y — u a Y, a — Y, ct . Einstein's field equations for (0) and (||) reduce to 

Y 2 = ™ + E, (3) 

4ttG _ M' 
c 4 P ~~ y 2 Y' ' 



where 



47rG , , 

M = — m (r), (5) 



and m(r) is an arbitrary function (with units of energy) that emerges as an "integration constant" . The form of the 
integrals (ie antiderivatives) of (|^) depends on the sign of the function E(r). Since (|3|) is the dynamical equation for 
LTB solutions, the usual convention classifies their possible dynamical evolution according to the three possible signs 
of E, as "parabolic" (E = 0), "elliptic" (E < 0) and "hyperbolic" (E > 0). Two possible situations arise: either 
E has a specific (positive, negative or zero) sign for the full domain of regularity of r, or only for an open subset, 
r\ < r < r2, possibly changing sign for other subsets of this domain. Whenever is necessary we will distinguish these 
two situations by using the terms "parabolic, elliptic or hyperbolic solutions" (full domain) and "parabolic, elliptic or 
hyperbolic regions" (subset). A convenient terminology that includes genericaly both situations is that of "parabolic, 
elliptic or hyperbolic dynamics". 

Only for parabolic dynamics we can obtain a closed integral of (||) expressible in the form Y = Y(t, r). For elliptic 
and hyperbolic dynamics we have either an implicit canonical solution of the form ct — ct bb {r) = F(Y,r), were F is 
the integral quadrature of (H) and t bb is an arbitrary function (ie an "integration constant" ) , or a parametric solution 
of the form [Y(rj, r), ct(rj, r) — ct bb (r)]. The function t bb marks the proper time value corresponding to Y = for each 
comoving observer, hence it is customarily refered to as the "bang time" . The solutions of @ are 

• Parabolic solution, E = 0, Y > 0: 

Y 3 / 2 = ±^V2Mc[t-t bi (r)], i>U+), *<*»(-) (6) 

• Elliptic solution, E < 0, < Y < 2M/\E\: 



, s JY(2M- \E\Y) M ( \E\Y 

ct(Y,r) = ct bb - V V - , 1 1 ; + — — arccos 1 - 



E\ |S| 3 / 2 V M 

irM 



Y > 0, ct.. < ct < 



El 3 / 2 ' 



ct(Y,r) = ct bb + 



y/Y(2M -\E\Y) M 



\E\ \E\V 2 



2tt — arccos 1 



\E\Y 



M 

ttM 2-kM 



3 



M 

Y(v,r) = j- [1 - cosr,}, (8) 



ct(r),r) = ct bb + te^H I V ~ sin *> ] ' ( 9 ) 



|£;| 3 / 2 



• Hyperbolic solution, E 1 > 0: 



±c(t-tj = V 1 f ^ - arccosh \l + — ) , (ID) 

M 

Y(r,,r) = — [cosh v - 1 ], (11) 
M 

ci(?7,r) = ct 6b ± --^ [sinhrj - ?? ] , (12) 

where the signs ± in the parabolic and hyperbolic cases distinguish between "expanding" or "collapsing" dust layers, 
that is Y increasing (+) or decreasing (— ) as t increases, so that for any given layer r = const., we have t > t bb 
(expansion) and t < t bb (collapse). Notice that the sign ± is not necessary in the elliptic case, since dust layers bounce 
(ie Y — 0) as Y — 2M/\E\ and so both expansion and collapse are described by the same equations (^) and (||) and 
by each branch in the canonical solution (ffl). 

As we mentioned before, t bb (r) marks the proper time associated with Y = for each comoving layer (ie the "bang 
time" ) . The conventional interpretation for the other two free functions M and E is based on the analogue between 
the dynamical equation ([|) and an energy equation in newtonian hydrodynamics (see jjj), in which M and E are 
usualy understood, respectively, as the "efective gravitational mass" and the "local energy per unit mass" associated 
with a given comoving layer. An alternative interpretation for E has been suggested as the local "embedding" angle 
of hypersurfaces of constant t ( |Q, p]|, pi]]). The specification or prescription of M and E is often made in the 
literature by selecting "convenient" mathematical ansatzes that are somehow intuitively based on these interpretations, 
loosely based on newtonian anlogues. However, configurations without centers p3| have no newtonian analogue and 
so the interpretation of M and E in this case is still an open question. Also, there is no newtonian equivalent for t bh 
and no simple and intuitive way in which this function should be prescribed, though the gradient of this function has 
been shown to relate to growing/decreasing modes of dust perturbations in a FLRW background @ @ @ @, 
and thus the functional form of t bb can be suggested from this feature. 

III. NEW VARIABLES 

A useful formulation of LTB solutions follows by redefining the free parameters M, E, t bb in terms of an alternative 
set of new variables defined at a given arbitrary and fuly regular Cauchy hypersurface, %, marked by t = ti (all such 
hypersurfaces are Cauchy hypersur face s in LTB solutions M). In order to relate the free function E to the scalar 



3-curvature of %, we remark from (A3) that the Ricci scalar of these hypersurfaces is given by 



a relation expressible in terms of invariant scalars (see (A3) in the Appendix). We introduce 



K ee —E, (14) 

so that the sign of ^TZ in ( Jl3| ) coincides with the sign of K. From here onwards we shall systematicaly replace E by 
—K, so that the type of dynamical evolution will now be denoted by the sign of K as: parabolic (K = 0), hyperbolic 
(K < 0) and elliptic (K > 0). Bearing in mind (§, @ and Q, we consider the following functions: 

Y(r) ee Y(U,r), (15) 
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AttG , . 
—TPi(r) 



AttG , . 



M' 
Y 2 Y! 



(16) 



Y- Y! 



(17) 



as initial value functions denned in %, The subindex j will indicate henceforth evaluation along t = i,-. For the 
moment we will assume this hypersurface to be fuly regular, we discuss the corresponding regularity conditions in 
sections VA and VII. 



In order to recast (|3|) and its solutions in terms of the new variables it is convenient to scale Y with respect to Yi . 
By defining the adimcnsional scale variable 



V 



Y 

y: 



and using @, @, @ and @, equations (Q), @ and @ become 

- r 2 (y ^2 



ds^ = -c A dt l 



■dr 2 + Y 2 (d9 2 + 



sm 



(18) 



(19) 



2pL 



Y- Y! 

■I _l 

Y 2 Y> 



K. 



where 



Pi 

y 3 r' 



y'/v 

YllYi 



(20) 
(21) 

(22) 



M 



(23) 



K 

Y?' 



(24) 



The solutions for the evolution equation (^p|) arc those of (||), that is, equations @ to (|lj) expressed in terms of the 
new variables y, /i, k by means of (|l^), (|23|), an d by eliminating t bb by setting t = ti and Y = Yi into equations 
@ to @. The solutions are 



Parabolic solution 



1 ± -v^cft - *j ) 



2/3 



(25) 



• Elliptic solution 



ct(y,r) = cU + 



[2^-n] 1 ' 2 -[y{2 l i-n V )f 2 



,3/2 



arccos 1 



arccos [1 

i 1 



y > 0, ct bb < cU < 



7T \JL 
73/2 



(26) 



Cf,; 



2 f i-K} 1/2 + [y(2[i- K y)} 1/2 /x 



f 3/2 



arccos 1 



ny 



arccos (1 I — 2tt 
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7T Lb 2n LI 



y{r],r) = - [ 1 - cos 77 ] 

K 



ct(r),r) = cU + [(r) -sin 77) - (rji - sin 77,) ] . 



(27) 
(28) 
(29) 



i] = arccos 1 — 



rji — arccos 1 — 



(30) 



• Hyperbolic solution 



,1/2 



ct (y,r) = cti± [y^+\«\y)] 1/2 -l^ + \«\ ]-'- ± I 



3/2 



arccosh 1 



my 



arccosh 1 



(31) 



y{v,r) = [ cosh t) - 1 ] , 



(32) 



13/2 



(sinh 77 — 77) — (sinh % - rji) } , 



(33) 



arccosh 1 



\ K \y 



arccosh ( 1 H 



(34) 



where, as before, the signs ± distinguish between expanding layers (— sign, t bb in the past of ti) an d col lapsing ones 
(+ sign, t bb in the future of ti). Explicit forms for t bb in terms of /1, k, rji and U are given in section VII . 

It is important to remark that y{ti, r) = yi = 1, by definition, and so y = 1 corresponds to t = ti. This is clear for 
the parabolic ( p5|) and hyperbolic solutions ( |3l|) and (33). However, it is important to notice that cti < ttll/k 3 ^ 2 in 
( p6| ) and so the initial hypersurface must be defined in the expanding phase in the elliptic solution. In this solution 
the correspondence between t = ti and y = 1 only holds in the expanding phase of (^6|) but not in its collapsing 
phase. The explanation is simple: in general the expanding and collapsing phases of the elliptic solution are not 
time-symmetric, as in the homogeneous FLRW case, therefore y — 1 will also be reached in the collapsing phase but 
it will not (in general) correspond to a single hypersurface t = constant. Figure [I] illustrates this feature. 



IV. VOLUME AVERAGES AND CONTRAST FUNCTIONS ALONG THE INITIAL HYPERSURFACE. 



The definitions (16) and ( |i7| ) of the initial value functions pi and ^IZi lead to the following integrals 

AttG 



M(r) 



(35) 



K{r) = 



2Y r 



(z) 1liY?Y' dr 



(36) 



where the lower integration limits in these integrals have been fixed by the conditions M(r c ) = K(r c ) = 0, where 
ri = r c marks a Symmetry Center^ (to be denoted henceforth as "SC" or as "center" , see section |V A and the 



1 See section V A for more details. See the Appendix for a formal defintition of a Symmetry Center 
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Appendix). For the remaining of this paper we shall assume, unless explicitly specified otherwise, that the spacetime 
manifold admits at least one SC and so the integrals (|35| ) and ([56]) become functions that depend only on the upper 
integration limit r. 

Since 4ir J Y^Y- dr is the volume associated with the orbits of SO(3) in the hypersurface % (though it is not a proper 
volume related to (A12)), equations (|35|) and (B6h suggest then the introduction of the following volume averages 



(PiKr) = 



( (3) ^)(r 



f ptYfYldr _ 3m(r) 
f r WlZiY?Y/dr _ 6K(r) 



fY?Y/dr 

In terms of these volume averages the variables \x and k in ( p3| ) and 

AirG 



n 7 

I) become 



P 



3c 4 



{Pi), 



i( (3) ^>, 



(37) 
(38) 

(39) 
(40) 



Setting pi and ^IZi to constants in equations ( pT\ j and ( p8[ ) yields = (pi) and ^ 3 ^i = (^IZi), therefore the 
comparison between pi, ^IZi and (pi), (^ s 'TZi) must provide a "gauge" that measures the inhomogeneity of pi and 
^TZi in the closed interval r > r c . This motivates defining the following "inhomogeneity gauges" or contrast functions 



Am) _ 



(pi 



- 1, 



A 



(fc) 



A 



(41) 
(42) 



These functions can be expressed in terms of M, K, p and k and their radial gradients. For example: 



• A^ and A^ in terms of fi, M, k, K 



A m ) 



A 



47rG Pi 
3c 4 # 

,/.) = { ^n± 

6 K 



4ttG Pt Yf 
3c 4 M : 



6X : 



(43) 



• A< m) and A\ K ' in terms of fi', M\ k', K' and Y- 



A k ) 



. ( m ) 



1 + 2 A 



(fc) 



M'/M 
K'/K 



p'/p 

K 1 1 K 



= A 



(m) 



2Y//Y;' 3Y//Y 



(44) 
(45) 



Another useful relation is 



M' 



The interpretation of A^™ 1 ' and A\^' as contrast functions will be discussed in section VI. Using ([I4[), ( |45[ ) and ( |46| ) 

we can eliminate all radial gradients like //, k' (or M', K') in terms o f A,- "- 1 and A J - fc ' ) . Since these gradients appear 
in the evaluation of T from ( |22| ) and the solutions of (20) (see section VII), it is certainly useful to be able to relate 

an important quantity like T to initial conditions expressed in terms of quantities like A^" 1 "* and A^ , quantities that 
have a more intuitive and appealing interpretation than M' and K' . 



3 K' 
2 ~K 

A k ) 



p_ 

P 



3 k' 

2 ~K 



Yi 



^(™) _ 3 ^(fe) 



(46) 
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V. GENERIC PROPERIES OF THE NEW VARIABLES 



We discuss in this section the regularity properties and generic features of the new variables, as well as their relation 
with M, K defined by (f{7|), (|38|). The averages and contrast functions, ( f4l| ) and (ff2|), are examined in section VI. 
Frequent reference will be made to expressions presented and summarized in the Appendix. 



A. Initial hypersurfaces. 

A hypersurface % is the 3-dimensional submanifold obtained by restricting an LTB spacetime to t = ti. This 
submanifold can be invariantly characterized as the set of rest frames of comoving observers associated with proper 
time t = ti (see the Appendix). Its induced 3-dimensional metric is 



ds 



Yf- 



1 — k Y/ 



dr 2 + Y? [dO 2 + sin 2 6»d(/> 2 ] 



(47) 



where < 9 < n and < < 27r and r extends along its domain of regularity (see subsection V C and the parragraph 
below). For a regular %, all invariant scalars like ( CO), (A4) and (A5) resticted to t = ti, or invariant scalars computed 
directly from (flTn, must be bounded. The mctric(M) is also subjected to the regularity condition (A7) which now 
reads 



K 



kY > 



l(3) 



&1Zj 

(^iz % 



K < 1 + A 



(fc) 



(48) 



where the equality must coincide with a zero of Y( (see below) and we have used (14), (|24|), (p8|), ( [40|) and (|42|). The 
hypersurface % can be foliated by the orbits of SO (3) marked bt r = c onst ., thus the regular domain of r is then 
limited, cither by a SC, or by those maximal values r = r max defined by (All) restricted to t = ti 



r Yjdr 



YJ dr 



1/2 



Y'dr 



1 + A 



(fc) 



(3) 



KiY? 



1/2 



oo as 



(49) 



where JC and k have been eliminated in terms of ^TZi and A^ by means of ([l0|) and (|4|). These maximal values could 
correspond to r — > oo or r — > ±oo, but this is not necessary, as it is a lway s possible to define the radial coordinate so 
that maximal values r max correspond to finite values of r (see section VC). 



• Regularity and differentiability. 

From their definition in (|l6|), ( [l7|) and (|l5|), the functions pi, ^TZi, Yi satisfy the same regularity and differ- 
entiability conditions as p^^lZ, Y in the radial direction (see the Appendix): Yi > must be continuous and 
Pi > o, ^n, can be piecewise continuous but with a countable number of finite jump discontinuities. The 
average functions, {pi), y^'lZi), constructed from intergration of pi, ^TZi, will be, in general continuous (like 

M and K). The contrast functions, A^ and A^ k \ explicitly contain pi and ^IZi (from (43)), hence they are, 
at least, piecewise continuous. 

• Zeroes of K and Y' . 



The zeroes of Yi and Y( are also zeroes of Y and Y', though the converse is not (in general) true. Because of 
( |l5| ) all regular zeroes of Y and Y' , associated with a comoving r (symmetry centers and turning values of Y') 
will be shared by Yi, but a regular % must avoid the zeroes of Y and Y' that do not correspond to a comoving 
r (singular zeroes). Therefore, we have along % 



— Symmetry Centers, SC.. 

Let r = r c mark a SC (see the Appendix). The functions Yi, pi 
and ([lj]), must satisfy 



(3) 



TZi, from their definitions in (|T5|), ( |l6| ) 
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Yi(r c )=0, Y i '(r c ) = S i ?0, *i«$(r-r c ) 
p t {r c )>0 p'i(r c )=0 p t ^ Pl {r c ) + lp'l{r c ){r-r c ) 2 

(50) 

WKi(r c ) yt ^^(r c ) = « ^TZ t (r c ) + ± (3) 7^'(r c )(r - r c ) 2 

(3) ^(r c ) = (3) ^-(r c ) = (3) ^i(r) = 0, parabolic dynamics 

where Si — S(ti) = Y'(ti, r c ) > and the fourth subcase above denotes the case of parabolic dynamics in 
which the condition K = k — holds in an interval r c < r, thus implying that ^TZi(r) = for all the 
interval. The central behavior of M, K in (|3^) and (|36|) is given by the leading terms 



K « ~ ^TZ^Sfir-r^, (51) 
o 

From (|23), (^3), (39), (f40|), (41) and (^2|), the value of /x, k and the remaining initial value functions at 
r = r c are given by 

K r e) = ^rPr^c), n(r c ) = i (3) ft 4 (r c ), 
(Pi)(r c ) = Pi(r c ), ( (3) ^>(r c ) = < 3 >7^(r c ), A^(r c ) = A 2 (fe) (r c ) = 0, (52) 

Signs and turning values of Y( . 



From (15) and comparing (|4|) with (16) and ( ]l3| ) with (|l7|), it is evident that regularity of p and P i requires 



the sign of Y( to be the same as the sign of Y' . Hence condition (A6) can be restated as 

Pi > «*■ sign(M') = sign(F/), (53) 

so that regularity of pi, ^IZi and of the metric (^7|) leads to the condition that all the following equations 

Y[{r*) = 0, M\r*) = 0, 1 - K{r*) = 0, (KY)' = 0, (54) 

must have common zeroes (turning values) , r = r* , of the same order in r — r* . Failure to comply with 
(pi|). not only reflects "bad" choices of pi, ^IZi, but leads to surface layers in r — r* associated with a 
discontinuous extrinsic curvature of the submanifold r = r* , [0 , |ll| , jL7| , JlJ] . 

B. Comparison of y vs Y and r vs Y'. 

Since the zeroes of {Y = 0, Y' = 0} do not (in general) imply {K; = 0, Y! — 0}, but all comoving zeroes of 
{Y = 0, Y' = 0} are also zeroes of {Yi = 0, Y( = 0}, those zeroes of {Y = 0, V' = 0} that are not common to 
{Yi = 0, Y( = 0} must correspond to curvature singularities (central or shell crossing). This removes the ambigueity 
of identifying these singularities by the constraints Y = and Y' = 0, placing this identification on y and T. From 
its definition in (|l8|), it is evident that y(ti,r) = 1, while its central behavior is given by 

V(t,r c ) = y(t,r) « ^ [ 1 + 0(r - r c ) + 0(t - U)} , (55) 

where S(t) — Y'(t, r c ), so that (in general) y does not vanish at a SC. This means that y = implies Y = 0, but the 
converse is false, since F = does not imply y = (for example at r = r c ). Therefore, the coordinate locus for a 
central singularity (either big bang or big crunch) becomes 

y(t, r) = 0, Central singularity. (56) 

an unambiguous characterization, as opposed to Y — which can be associated either to such a singularity or to a 
SC. In fact, as mentioned in the previous section and from the form of the metric (|l^), the function y plays the role 
of a local scale factor, in contrast to Y defined non-localy in terms of the proper surface of the orbits of SO (3) around 
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a SC that might not even be contained in the spacetime manifold. The difference between y and Y is illustrated in 
figures |^ displaying parametric 3d plots Y(t,r) and y(t,r) for various LTB models. These figures clearly show that 
y(ti) — 1 for all layers. 

From (AS), the gradient Y'/Y is singular at a SC, however, from ©, (H|) and (H), we have 



r(t, r c 



1. 



■S ; 



S3(t)- 



0. 



3S(f) 
S(t) 



while the coordinate locus for shell crossing singularities can be given as 

T = 0, pi > 0, Shell crossing singularity. 



(57) 



(58) 



a characterization that avoids the ambigueity of Y' = 0, since, as mentioned in the Appendix, Y' might vanish and 
change sign regularly for a single r = r* ^ r c if (A9) holds, preventing surface layers. This situation is now taken 
care by conditions d5~3) acting on Y!. 



C. The function Yi m . choice of topology and radial coordinate. 

From the common geometric meaning of Y and Yi (see Appendix), it is evident that the homeomorphic class 
( "open" or "closed" topology) of the regul ar hy persurfaces of constant t follows from the homeomorphic class of Ti . 
The latter can be determined by applying flAl2Q to ( f47f ) for a given choice of Yi (r) , along the domain of regularity of 
r limited between SC's (real roots of Yi(r) — 0) and/or a maximal value r max defined by (|49|). Because we can always 
rescale r, it is always possible to select this coordinate so that Yi takes the simplest form, hence the freedom to select 
Yi can also be the freedom to find a convenient radial coordinate for a desired LTB configuration. As mentioned in 
the Appendix, the compatibility of any one of the possible choices of homeomorphic class with any one of the three 
choices of dynamical evolution (parabolic, hyperbolic or elliptic) is only restricted by the regularity conditions dig), 
(|54|) and (|i^). We ellaborate these issues below looking separately the case with one and two SC's. 

• One symmetry center. 

Since we can always mark the SC by r — 0, then 3^(0) = and if Yi is a monotonously increasing (one-to-one) 
function (X/ > 0) for all r, we can then choose the radial coordinate so that Yi — Sir, where Si is a characteristic 
constant length. This is equivalent to using Yi/ Si as a new radial adimensional coordinate, but we can also set 
Yi = Si f(r), for any monotonously increasing function with /(r c ) = 0. The range of r is < r < r max , so that 
proper distances £i increase as r grows and £j — > oo as r — > ?' max . A useful coordinate choice is / = tan r, so 
that maximal values associated with Yi — > oo can be mapped to r = ir/2. Fulfilment of conditions (E^, ( |5~i| ) 
and ( |l9|) is trivial for parabolic and hyperbolic dynamics, but for elliptic dynamics we must have K — > 1 as 
r — > r max . The homeomorphic class is R 3 . 

• Two symmetry centers. 

Let the centers be marked by r = r Cl and r — r C2 . Then Yi{r Cl ) = Yi(r C2 ) = 0, so that Y[ vanishes for some r* in 
the open interval r Cl < r* < r C2 . In this case, it is not possible to rescale Yi with r because Yi is not a one-to-one 
function of r. However, depending on the model under consideration, it is always possible to select the radial 
coordinate so that Yi = Si /(r), where f(r) is a C° function satisfying f(r Cl ) = f(r C2 ) = and f'{r*) — 0, so 
that pi ', ( - 3 - > IZ'i, M', K' ', 1 — K = 1 — nYf all vanish at r* and have zeroes of the same order on r — r*. Since 
the hypersurface t = ti has a finite proper volume, £i evaluated from one SC to the other (the maximal range 
of r) must be finite. Although it is possible to set K < (parabolic and hyperbolic evolution), only for elliptic 
evolution it is possible to fulfil condition I p4\ ) in order to avoid surface layers in r = r* (see |l6| and p8|). 
The hypersurface % is homeomorphic to a 3-sphere, S 3 . 



D. Invariant expressions. 



The volume averages and contrast functions introduced in section |y| were constructed by averaging pi and ^IZi 
with the volume (4/3)TrY i 3 = in fY?Y-dr, a quantity that can be related to the proper area, 47rF i 2 , of the group 
orbits of SO (3). However, it can be argued that these volume averages and contrast functions are coordinate dependent 
expressions. We can show that these quantiti es a re equ ival ent to coordinate independent expressions given in terms 
of invariant scalars evaluated along %. From (A4) and ( A5) applied to t = ti, together with (f4l|), we obtain 



10 



iirG , . 8ttG r , r . 

— r" \Pi) = —j- Pi + 6 W2 \i = TZi + o[ip2\i 



1 + A (m) = 4nGpi = Ul (W\ 

where TZi and [^2]i ar e the 4-dimensional Ricci scalar and t he conformal invariant 4*2 evaluated at t = tj. The 
evolution equation © applied to i = ij, together with (|aJ), ([A^) , ©, ©), © and (f§) yields 



1 + AW = ggj + llQ^lgjg] = C 2 (3)^ 

C 2 % + 6c 2 W,-((T ! + ie ! ) 2 c 2 ( 3 )^ - 12o-i(a-i + |e0' 



where cr^ and Oi are the kinematic invariants appearing in (Al) and JA2J ) (the expansion scalar and the degenerate 
eigenvalue of the shear tensor) evaluated along %. The expressions (|59|) and (60) are not the definitive "covariant 
interpretations" for the average and contrast functions. They have been presented in order to emphasize that these 
quantities can be expressed in terms of invariant scalars and so can be computed in any coordinate system. An 



alternative volume averaging could have been attempted using the proper volume (A12), leading to quantities like 



5 Pi {l-K)-V*Y?Y!dr fWKijl-KrWYfY'dr 
J(l-K)-^Y?Y>dr ' f(l-K)-V*Y?Yfdr ' 1 ' 

instead of the volume averages in (|37| ) and (|38|). Notice that only for parabolic dynamics all these averages coincide. 
We have chosen to use (|37]) and (|38|) instead of (|6l]) mostly because the former are mathematicaly simpler than the 
latter and because they follow naturaly from the field and evolution equations. For example, the functions M, E = —K 
or /i, k that appear in the solutions of (^) or (^0|) immediately relate to (|37]) and (^), but their relation with (^H]) is 
(except for parabolic dynamics) mathematicaly difficult, making it very hard to accomodate quantities like ( |6 1[ ) to 
the dynamical study of the models. Since all the invariant scalars defined in by (^) and (|6^) can be generalized 
for arbitrary times t ^ ti, the volume averages and contrast functions can also be generalized in this manner. For 
example, the fact that M = M(r) and K = K(r), and so these functions are constants for every comoving oberver, 
follows (from ( p7| ) and (p8|)) from the conservation (along the 4- velocity flow) of quantities like (p)Y 3 and (^TZ)Y 2 , 
were (p) and (™7?.) arc the generalization to arbitrary times of (pi) and (^TZi). This and other conservation laws 
are being examined elsewhere p9| . 

E. FLRW and Schwarzschild limits, matchings and hybrid models. 

It is a well known fact that LTB solutions contain as particular cases FLRW cosmologies with dust source (ho- 
mogeneous and isotropic limit) and the Schwarzschild-Kruskal spacetime (vacuum limit). We examine these limits 
below 

• FLRW limit . If for all the domain of regularity of r we have p z = pf } and ^TZ % =^ Tzf \ where pf ] > and 
( 3 )^ 0) are constants, then equations @, Q, @ and @ imply 



Pi = M = Pf\ =► aW=o, M = *£ P f>Y? 

W-Ri = (^TZt) = ( 3 >ftj°>, Af } =0, K = 1(3)^^, (62) 

Since a and k in (123) and (E4J) are now constants, y (determined by (EQ)) must be a function of t only, and so 



from (18) and (|22|j7we have Y = 1 and Y = (S{t) / Si)Yi, identifying y = S(t)/Si with an adimensional FLRW 

scale factor. Rest mass density and scalar 3-curvature take the FLRW forms: p — pf^ (Si/ S) 3 and ^TZ =^ 

TZf ] (SJS) 2 , while the metric © becomes a FLRW metric if we identify Y t = S t r so that (1/6)( 3 )^ 0) 5| 

becomes an adimensional curvature index. Equations ([32|) reafirm the role of inhomogeneity gauges for A^ 

and A\ k \ making also evident that initial conditions with homogeneous pi and {3) TZi leads only to FLRW 
evolution. 
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• Schwarzschild limit . For an intial density given by a normalized Dirac delta (a point mass) 

(H = P? ) 6(Y?) = pV> s!s(r3) 



(63) 



we obtain: m = p' ' Sf/3 from ([35]), so that we can identify the "Schwarzschild mass" with 2M given by (||). 
The metric (|l|), together with (|^)with M = const, describes the Schwarzchild-Kruskal spacetime in terms of 
the worldlines of geodesic test observers with 4- velocity u a . The function ^TZ is the scalar 3-curvature of the 
hypersurfaces of simmultaneity of these congruences. 

If only one of pi or ^IZi is constant (but nonzero) in a domain of r containing a SC, then only A^ m ' or A^ fc ' vanishes 
in this domain. For hyperbolic or elliptic dynamics, these are special cases of initial conditions but do not correspond 
to the FLRW limit (for parabolic dynamics ^TZi = and so the FLRW limit follows from pi — pf- 1 > 0). The cases 
in which pi and ( 3 '7Zi are, in general, not constant but become constant (even zero) only in an interval of the domain 
of r can be also understood as initial value functions of "hybrid" spacetimes obtained by matching a given LTB region 
(with pi, ( 3%> lZi not constants) with a section of FLRW or Schwarzschild spacetimes in another part. These matchings 
can also be defined among various types of LTB (non-FLRW or non-Schwarzschild) regions (see jll| , JlJ] ) . We will 
not examine these hybrid spacetimes in this paper (see [p9[). 



VI. LUMPS AND VOIDS IN THE INITIAL HYPERSURFACE. 
A. Initial density lumps and voids. 

The relation of pi with (pi) and is straightforward. Consider the following auxiliary equations obtained by 

applying integration by parts to ( p7| ) 

1 f r 

Pi - (pi) = ^73 / Pi Y i dr 

In all intervals of r containing a SC, p[ ^ implies A^™' ^ 0, except at the SC itself where both p\ and A'" 1 '' must 
vanish (from © and ©). We still have A| m) ^ 0, even if pi and/or p\ vanish in a subset of this domain that 



contains no SC. With the help of ( p7| ) and (|6J), we discuss below the cases with one and two symmetry centers. Unless 
specihcaly mentioned otherwise, we assume that either p\ > or p\ < throughout the interval under consideration. 
The discussion below is illustrated in figures 

• One symmetry center. See figures |a and |b. 
In this case a density "lump" or "void" is characterized by pi having a maximum/minimum at the center 

(m) 



worldline r — r c . Hence, for a lump or void, we have (from (64)) the following possible ranges of A^ 
initial lump: p\ < 0, p 4 < (p t ), => -1 < A| m) < 0, 

initial void: p\ > 0, p, > (p l ), => A| m) > 0, (65) 

where the lower bound on A^™"* for a lump is A^ = —1 if pi vanishes for r > r c (from(^l|)). 
Two symmetry centers. See figures |= and |d. 

From (|0|), (|l|) and (||), all of Yi, M, M', Pl \ Yi, pi - (pi) and A^ m) must vanish at both SC, r = r Cl and 
r = r C2 , with zeroes of the same order in r — r Cl and r — r C2 . This means that there must exist a turning value 
r = r* in the open domain r ci < r < r C2 for which Y(, M', pi ' all have a common zero of the same order at 
r = r* . Besides these conditions, pi and Yi must satisfy 

M(r C2 ) = ^ f " Pl Y?Yldr = (66) 
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Therefore, for every r — a in the open domain r ci < r < r C2 we have 

Pi Y?Y(dr = - f" PiY?Y(dr = [ Pi Y?Y?dr (67) 



Since Yi vanishes at both SC's, this means that, for whatever pi we choose, the average (pi) evaluated from one 
center (say r = r ci ) to r = a must be the same as that average evaluated from the other center (r = r C2 ). Also, 
following the usual interpretation of M as the "effective mass inside" a comoving layer r, this "effective mass" 
for the whole hypersurface (from r = r Cl to r = r C2 ) is zero. 

We have two possible regular configurations (see figures and [|d), depending whether pi is: (i) a local maximum 
at each center (each center is a density lump) and (ii) a local minimum at each center (each center is a density 

void). For each case, A l - m '* vanishes at the centers, while (from ([34])) and using the same arguments as in the 
case of one center, we have for the cases (a) and (b). 

an initial lump in each center: pi < (pi), => — 1 < A^ < 0, 

an initial void in each center: pi > (pi), => A^™ 1 ' > 0, (68) 

If allow pi ' to change its sign in a given domain containing a SC, we obtain an initial density that can be associated 
with "density ripples" (see figure ^a). In this case we have — 1 < A^ < for all r for which pi' is negative and 
Aj > for positive pi ' . 



B. Initial curvature lumps and voids. 

Unlike pi, the intial 3-curvature can be zero, positive, negative or change signs in any given domain of r. This 
makes the relation between ^TZi, (^TZi), K and A^ more complicated than in the case of pi. Also, the regularity 
conditions ([l8]) and (|54|), as well as the distinction between the type of dynamics (parabolic, hyperbolic or elliptic), 
all involve specific choices of ^IZi and K but not of pi and M. 

The following auxiliary relation, analogous to (|6^), follows from integration by parts of ( |38| ) 

W-Ri - (Wlli) = r^Tli'Y^dr, (69) 

M Jr e 

In all intervals of r containing a SC, ^TZ^ ^ implies A^' ^ 0, except at the SC itself where both quantities vanish 
(from ( po|) and (|5fj|)). If ^TZ{ = for a domain of r containing a SC, then (^IZi) = K = and this domain is a 
parabolic region (or a parabolic solution if it is the full domain of r). However, we mi ght have (3) ^ = and/or 
^TZi ' = in an open subset of the domain of r containing no SC's without K or A^ vanishing. We look at various 
cases below. The discussion below is illustrated by figures f|, || and |]b. 

• ^TZi > for a domain of r containing a SC. 

The relation between ^TZi, (^TZi), K, and A^ is analogous to that between pi, (pi), M and A^ (see figures 
^a, |]b and ^c). Since (^TZi) and K are everywhere positive, the dynamics of the dust layers is elliptic. 

The results for the cases of curvature "lumps" or "voids" (sign of ^IZ'i) and/or when there is one or two centers 
are analogous to those discussed earlier for pi : 

curvature lumps: (3) 7^ < ( (3) ^), -- < Af ' < 0, 

curvature voids: {3) TZ t > ( (3) ^), A[ k) > 0, (70) 

where the minimal value —2/3 follows from K > and from the regularity of the left hand side in (^). Notice 
that condition (|^) is not trivialy satisfied. We examine below the cases with one and two SC's. 
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One SC. (See figure |a). 

From (Q) and (pi]), the functions Yi and K vanish at the center and increase as r — ► r max , but (H) requires 
that K < 1 for all the domain of r. Condition (p4) can be avoided by selecting Yi = Sir so that Y( never 



vanishes. Necessary and sufficient condition for complying with (48) and ( [49|) are 



i r 

m 1Z, ^0 and K = — J ( 3 >KiY?Y-dr -> 1 as r -> r max , (71) 

so that ti — > oo as r — > r max . The model examined in Q is a particular case of this type of configuration. 
- Two SC's. (See figures |b and|c). 

In this case all the functions K, K', ^TZi — (^TZi) and A\ must vanish at both centers. Also, 
(3) 72-, ( (3) 7£i)', K' must all vanish for a value r* in the open domain r Cl < r < r C2 where Y( vanishes. In 
order to satisfy (^8|) and (|54|), the functions ^TZi and Y t to be selected so that the following conditions 
hold 

[ ° 2 WlZiYfYtdr = 2K{r C2 )Y t (r C2 ) = 0, 

Jr ci 

®1liY?Yfdr = K(r*) = 1, (72) 

The graphs of ^TZ i} {^1Zi),K and A\ k) are shown in figures [l| and qualitatively analogous to those of 
Pi, {pi), M and A^ in figures ||. The function if must vanish at the SC's (like M) and comply with J7^). 

< for a domain of r containing a SC. (See figures |a and |b). 

In this case we have (^IZi) and JC negative for all the domain of r, condition (^) holds trivialy and dust layers 
follow hyperbolic dynamics. However, the characterization of "lump" or "void" in terms of the sign of ' 3 '72^ is 
the inverse to that when ^TZi is nonegative: lump and void are now respectively characterized by ' 3 '7^ > 

and ^TZi < 0. Besides this point, the range of values of A l - fc ' ) for lumps and is the same as for the non- negative 
case in (70). 



^TZi changes sign but K and (^IZi) do not. (See figure |6Jt 

since (^1Zi) and K are obtained from an integral involving ( 3 '1Zi evaluated from a SC to a given r, it is quite 
possible to have situations in which ^TZi vanishes or pases from positive to negative (or viceverza) but (^TZj) 
and K remains positive (or negative) at least for an interval of the domain of r not containing a SC, so that 
the model remains elliptic or hyperbolic in this domain. The situation in which (^TZi) or K vanish for a given 
r = r* in an open domain of r that excludes a SC is qualitatively different and more complicated, since the 
dynamic evolution of dust layers (governed by (||) or (po|)) passes from elliptic to parabolic or hyperbolic (or 
viceverza). This case will be examined in a separate paper ]29|. 



VII. THE FUNCTION F, SHELL CROSSING SINGULARITIES AND BANG TIMES. 



The function T defin ed i n ( 22|) is ess ent ial for calculating relevant quantities characterizing LTB solutions, such as 
p, O, <t, ^72. and ip2 ((Al), ( A2[), (A2), ( |A5p and (^l|)). The exact functional form of T is also required in order to 
find the restrictions on initial conditions that fulfil the important regularity condition: 



T > for pi > and K ^ 1, 



(73) 



with r = only if pi = and/or 1 — K = 0. This condition garantees that: (a) the metric coefficient g rr is regular, 
(b) p > for a given choice of pi > and (c) provides the necessary and/or sufficient restrictions on initial value 
functions in order to avoid shell crossing singularities (these conditions are summarized in Table 1). We show in this 
section that the conditions for avoiding these singularities based on j7^ ) are equivalent to those found by Hellaby and 
Lake ^ (see also (l6) and 0). 

The function T can be obtained in terms of y and the initial value functions /i, k, A|™', A^ by eliminating y' /y 
from an implicit derivation with respect to r of equations riga), (p6|), (p7|) and (|3l|). For hyperbolic and elliptic 
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dynamics, T in terms of r\ follows from implicit derivation of (|29j) and (|33|), using Q28h, (|32|), ( |30| ) and ( |34| ) in order to 
eliminate 77' in terms of y' /y. We provide below the most compact forms of this function for parabolic solutions (in 
terms of y) and for hyperbolic and elliptic solutions (in terms of 77) 



T = 1 + A| m) - parabolic, (74) 



A 



(m) 



r = 1 + 3(AH-Af 3 ) 



with P(rj) and Q(r)) given by 
• Hyperbolic dynamics 



1 - 



P 
P 



- 3P[Q-Qi] [ A[ m) - I Af 5 ) , hyperbolic and elliptic, (75) 



P(V) = 



sinh 77 \J 2 + x y 



(cosh 77 — 1 ) 2 (iet/) 3 / 2 ' 



Q{i]) = sinh 77 — rj = \Jxy(2 + xy) — arccosh (1 + x y), (76) 



V2 



Pi = P(Vi) = x3/2 1 Qi = Q(Vi) = v^T^Tx) - arccosh (1 + x ), 

\k\ 

?/ = arccosh ( 1 + x y ) > 0, Vi = arccosh ( 1 + x ) > 0, a; = — , (77) 



• Elliptic dynamics 



p(ri) = sin V = Vxy(2 - xy) = ± V2 - a;y 
( 1 — cos rj) 2 (xy) 2 (xy) 3 / 2 ' 



Q(v) = V — sm V = arcsin ( 1 — x y ) — y/xy(2 — xy), (78) 



V2 



f» = P(Vi) = V ^3/2 ~ > Qi = ( 3( ? 7i) = arccos(l-a;) - \Jx(2- x), 



77 = arccos( 1 — xy), < rj < 2tt 
r\i = arccos ( 1 — x ) , < r\i < n, x 



K 



(79) 



where all functions of 77 and r)i have been constructed with Q30|) and 
were eliminated in terms of A| m ^ and by means of equations (| 



) and the gradients of the form \j! j \x and k'/k 
(|f) and ©. Notice that the ± sign in @ 
follows from the fact that sin 77 becomes negative in the collapsing phase (n < 77 < 27r), whereas sinh (77) is positive 
for all 77 < 0. It is possible to use ([77|) and ([79|) in order to transform T in (75) as a function of y instead of 77, but 
the obtained expressions are more complicated than (f75|). Also, the form of T = T(y,r) for elliptic solutions consists 
of two branches, making it very impractical for most purposes. 

Since the central and shell crossing singularities are (in general) not simmultaneous with respect to t, the existence 
of regular initial hypersurfaces t — t{ is not a trivial issue. Such hypersurface must avoid the coordinate locii of shell 
crossing and central singularities. The conditions for avoiding shell crossings will be obtained below, not just for t = U, 
but for all the time evolution, but assuming that the no-shell-crossing conditions have been met, the hypersurface % 
must still avoid the central singularity. The conditions for this can be given straight away as 



U 



tbb > 0' f° r expanding layers, 
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U — t bb < 0, for collapsing layers. 



(80) 



These regularity considerations provide an insight on the relation of the "bang time" t bb with the initial value functions 
Pi, ^TZi and their corresponding average and contrast functions. From the forms of t bb that will be obtained in this 
section and displayed in figures [?], the fact that t' bb ^ leads to an r dependent "age difference" for the comoving 

observers, a difference whose maximal value is very sensitive to A^ and A^ k \ and so to the ratios of maximal to 
minimal pi and ^TZi (but not to local values pi or ^TZi). This feature is illustrated in figures 0. 

The fulfilment of conditions (|7^), as well as a comparison with no-shell-crossing conditions of Hellaby and Lake 
|^5| , is examined below for the cases with parabolic, hyperbolic and elliptic dynamics. 



A. Parabolic solutions. 

From (f74|) and (|7|), shell crossing singularities are avoided if 



A.( m ) 



(1+aW)^_ A W>0 y^ > , (81) 



1 + A< m) ' 



1 < A< m) < 0, Pi > 0, (82) 



Therefore, necessary and sufficient conditions for avoiding shell crossings are 

-1 < a| to) 

conditions implying that we have initial density lumps. 

Notice that condition ( [82] ) will not hold, in general, for density voids (A^ > 0). However, since A| m) /(l+A t H)< 
1 and given a choice of a regular t = i,-, condition (^2|) is satisfied for initial lumps or voids and an evolution free 
of shell crossings is always possible for all y > 1 (or t > ti in expanding configurations). This fact allows for a well 
posed intial value study of LTB parabolic solutions with density voids, a feature that has been exploited in papers ( 
||, Q) dealing with LTB solutions in which shell crossings do occur, but their coordinate locus can be confined to 
earlier times t bb < t < U (see figure |TT|a) . 

The no-shell-crossing conditions of Hellaby and Lake ( p5[ , |l6| ], [p"7[) for the parabolic case and for expanding 
layers are given by 

± Y' > =*> {±M' >0, ±ct' bb < 0}, (83) 

where the ± sign accounts for the possibility that Y 1 might change sign regularly (only if M' vanishes). However, the 
equality of the signs of Y' and M' in (B3) follows from (p3|) and and m! = pY 2 Y' = piYfY[, and is already contained 



in (73). Regarding the sign of t bb we have, from setting y = in (|25|), the following expression for t 



2 , s 

ct bb = cti T ^-js=, (84) 



where the =p sign accounts for expanding (— ) or collapsing (+) layers. In order to compare with (|83]), we consider t' 
the case of expanding layers 

ct' = = -L ^A< m) , (85) 

3^/271 /it 7271 1 ' v ; 

This equation clearly shows that t' < for Y/ > only occurs if A^ m) < 0. Thus, the no-shell -crossing conditions 



of Hellaby-Lake (83) are equivalent to (p2|). However, we argue that phrasing these conditions in terms of pi, Yi and 

A„- m ' ) is more intutive than to do it in terms of Y', M' and t' . The function T for a regular model with parabolic 

^ I — — 1 - 1 

dynamics is qualitatively analogous to the plot displayed in figure |10|a. 

Given a choice of t = ti, it is evident from (B2h that the intial hypersurface avoids shell crossings, either for lumps 



or voids (as long as A^ is finite). Regarding the central singularity, the form of t bb in (|34|) implies that conditions 
( |so| ) are satisfied as long as p does not diverge, a condition that is, in turn, satisfied as long as pi and p\ are finite. 
As shown by figure [|, selecting unphysical initial conditions associated with "bad" choices of pi (diverging at the SC 
or for r > r c in a void) leads to a singular % for which Jso| ) is violated. As shown by figures 0a and 0b (the parabolic 
case is analogous to the hyperbolic one) , for reasonable initial value functions it is always possible to garantee that a 
regular hypersurface t — t t exists. In general, a "good" choice of M and t bb , complying with (|83|), will yield a "good" 
form of pi and Yi complying with ( [82] ) and (|S0|), however, we feel that selecting a 'good' form for pi is easier and more 
intuitive that choosing adequate forms for t bb and M. 
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B. Hyperbolic solutions. 



In order to examine the fulfilment of @, we assume that pi > and ^TZi < and rewrite (f75|) as 



r = 1 



A A 



(m) 



P A 



(k) 



3P r 
~Pi 



A A 



(m) 



Bi A 



(86) 



A( V ) 



P(v)Q(v), B( v ) 



A, 



Mm 



B t 



(87) 



where P(jf) and Q(rj) are given by ([76j). The fulfilment of (|73j) obviously depends on the behavior of A, B and P 
in the full range r\ > 0. We have P > with P — » oo as r\ — > 0, while in the same limit A — » 2/3 and P — > 0. 
Therefore, the following are necessary conditions for rt73) 

AAi m) - P,:A„ (fc) < 0, 



AA 



(m) 



PA 



(fc) 



> 0, 



(89) 



However, A and P satisfy < A < 2/3 and -1/2 < P < for all 77 > 0, with A -> and P -> 1/2 as 7/ -> 00. The 
functions .A, and Pi are bounded by the same values. Further restrictions follow by inserting these limiting values into 
(jS8[) and (|8j|). This leads to the following necessary and sufficient conditions for an evolution free of shell crossings 



> 0, 



(3) 



< 0, 



-1 < A| m) < 0, 



2 

-- < 
3 ~ 



< 0, 



(90) 



Since ( |88| ) and (89) imply (|90|) and the latter also implies the former, conditions (|90j) are, indeed, necessary and 
suffici ent. Notice that, just as for parabolic solutions, shell crossings cannot be avoided (in general, see subsection 
VII D| below) for positive contrast functions (initial density or 3-curvature voids). The functions T and p for a regular 
model with hyperbolic dynamics are shown in figures |To| a and |To|b. 

The no-shell-crossing conditions of Hellaby and Lake ( |l5|] , |16|] , |L7| ) for the hyperbolic case and for expanding 
layers are given by 



± Y' > =*> {±M' >0, ±K' <0, ±ct' < 0}, 



(91) 



where, as in the parabolic solutions, the ± sign accounts for the possibility that Y' might change sign regularly 
and we have used (|lj) in order to translate a positive E (in the original LTB variables) into a negative K. As in 
the parabolic solutions, the condition that Y', Y( and M' must have the same sign follows from (g3|), pi > and 
ml = pY 2 Y' = piY 2 Y(. Opposite sign for Y' and K' follows from (45) and ^IZi < (bearing in mind that K < for 
hyperbolic solutions). Regarding the sign condition on t' , we compute t bb by setting 77 = and y = in equations 
(|l|) and (||), leading to 



13/2 



2 \ 1 ^ 2 arccosh (1 + x) 
1 + 

X I X 



(92) 



where the ± sign accounts for expanding (— ) and collapsing (+) layers, 
expression given in terms of rji, we obtain 



ct'„ = 4^ 



3^7 

Yi \n\ 3 / 2 Pi 



A, A 



(m) 



Bi A 



(k) 



Considering expanding layers and the 



(93) 



but, because of (pq), the factor in square brackets must be negative or zero for all r\i > 0, thus implying that t' and 
Y( must have opposite signs. Looking at the maximal and minimal values of A^ and P^, together wit h (H ) and ( |45| ) 
leads then to conditions ( |9C|). Therefore, Hellaby-Lake no-shell-crossing conditions are equivalent to (|90D. 

Assuming that (|8^) and (^9|) hold, a regular % must satisfy conditions ( pcf ) in order to avoid the central singularity. 
Since the term in square brackets in ( |92"| ) tends, respectively, to and 1 as x — * and x — * 00, sufficient conditions 
for ( |30"| ) to hold are that both p and \k\ > be finite everywhere. These conditions are satisfied if | (3) ^| > and 
Pi are finite. As with the parabolic case, conditions (p0[) are violated only with very unphysical initial conditions 
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that might also violate the no-shell-crossing conditions (see figures f?]a and ^|b), since diverging pi and |^72.j| can 
be associated with for voids with large values of A^ and A[ k ^ . Hence, a regular hypersurfaces % always exist for 
physicaly reasonable intial conditions and a 'good' selection of functions M, E = —K, t bb leads to 'good' forms of 
^TZi, pi, Yi- However, we argue again, that finding a 'good' choice of the latter is easier and more intuitive. 

Shell crossing singularities will always emerge for initial conditions with density and/or 3-curvature voids (see figure 
fj]b), however for reasonable init ial co nditions characterized by finite p.; > and ^IZi < 0, we must have finite p > 



and k < in all % (see section VIH for examples). Thus, unless we choose pathological initial value functions as in 

figure H quantities like Pi and Qi defined in ((77|) will be finite and non-zero, leading to T — > 1 — (3/2) A^ as rj — > oo 

asymptoticaly, therefore a regular evolution of intial voids is possible in the domain t > f, as long as A^ < 2/3 
(for expanding layers) . An LTB model with initial voids and hyperbolic dynamics following this type of evolution is 
depicted in figure [TT|a. 

C. Elliptic solutions. 

As with the hyperbolic case, we assume that pi > and ^TZi > and rewrite (|7^) as 

r = 1 + 3(A| m) - Af) - 3PQ (aH - ~a|*>) - ^ [A t A { ™ ] - B, Af } ] , (94) 



A, 



1 



Pi Qi 



B, 



1 



PiQi 



p = p(th), Qi 



Q(Vi), 



where P(ri) and Q{rj) are given by (78). Necessar y a nd sufficient conditions for avoiding shell crossings can be obtained 
from ( |73| ) and ( |94| ) . Following Hellaby and Lake ]l5| , we explore the sign of the functions P and PQ near the big bang 
(j) ~ 0) and near the big crunch (77 w 2ir). As r\ — > we have P — * +00 and PQ — > 2/3, while PQ — > 2nP — > —00 as 
rj — > 27r, therefore necessary conditions for T > are 



A A 



(m) 



Bi A 



(fc) 



< 0, 



(95) 



2tt P ( A, (m) - I A| k) J + Aj A, (m) - P, Af 5 = (2tt - Qi) P l ( A,< m) - ~ A.| fc) J + A, (m) - A.[ fc) > 0, (96) 

Since the range of 774 is < < n, we have 1/3 < Ai < 1, < Pi < 1 and < Qi < tt, a third necessary (but not 
sufficient) condition that follows from (|)5|) and ([)6|) is 

A (m) _ 3 A ( fc) ^ ^ (g7) 



Following the arguments of Hellaby and Lake |15j , we prove the sufficiency of (B3) and (|96|) by assuming that 



1 + 3 A 



(A^-Af) = 3a (,) (Af - §A«) 



A M _ 3 (fc) = _/3gi) /^. A (m) _ ^ACtA (Q8) 



2 1 P 

for unspecified functions a(r]i), /3(?Ji), so that 

A. a("') _ n. aW 
r = _ . i ' ' [/3(a-PQ) + P], 

Pi 

and the fulfilment of T > requires ( j95|) and /3 > — P/ (a — PQ). However, from the forms of P and Q in (|94|), this 
latter condition can only be satisfied in the whole evolution range < r\ < 2ir if the functions a, (3 satisfy a > 2/3 

and (3 > l/(27r). Inserting these inequalities into ( p8| ) we obtain ( p6f ) together with A^ > — 1. Conditions ( |95| ) and 
( |96| ) are then necessary and sufficient, but we still need to test the compatibility of these conditions with the signs 
of A^ m ' and A^ (initial lumps or voids). Clearly, (E)5|) is compatible with either A^™ 1 ' < and > (density 
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lump and 3-curvature void) or with A^ m ^ < and A^ KI < (density and 3-curvature lumps). However 



(fc) 



is only 



compatible with the second choice. Therefore, the full set of necessary and sufficient conditions for avoiding shell 
crossing singularities are (p7|), together with 



: A a (fc) < 0, 



-1 < aH < 0, 



(99) 



and the two inequalities (^5|) and (|96|) which can be combined into the single expression 



2ttP A 



Am) 



- ^A (fe) 
2 1 



. p Q ( A (m) _ 3 A ( fc ) 



- A 



- A 



> 0, 



(100) 



A graphical representation of this condition is provided in figures [|a and [)]b. 

The no-shell-crossing conditions of Hellaby and Lake ( jla] , |lq| , |L7| ) for the elliptic case are given by 



± Y' > => {±M' >0, ±ct' < 0, ±ct' > 



-2irM ( M' 3K 



K 3/2 



M 2 K 



(101) 



where the ± sign accounts for the possibility that Y' might change sign regularly and (|l4|) has been used. As in the 
parabolic and hyperbolic solutions, the condition that Y', Y( and M' must have the same sign follows from pi > 0, 
(PI), ( p3| ) and ml = pY 2 Y' = piY?Y(. The same sign for Y' and K' follows from K > 0, (Ha) and (pq). Regarding 
the sign condition on t' , we compute t bb by setting r\ — and ?/ = in equations (|2q) and (E9J), leading to 



K 3/2 



arccos ( 1 — x ] 



1/2' 



< k <2p, 

We differentiate the term depending on r\i with respect to r 



< x < 2 



3Y/ 



ct:. = -TT- 



A, A 



(m) 



Si A 



(fe) 



(102) 



(103) 



where A4, Bi and Pj are the same functions given in (p4[). Comparing i' above with fl9q), it is evident that the 
Hellaby-Lake conditions (101) involving t' are the same as (D5h and (09), and so can be reduced to (IT' 



(100) if 



one translates these conditions in terms of A 
Assuming that 



(to) A (fc) 



A- ; and ty. 



|9q ) - (10C) hold, a regular hypersurface % must satisfy both conditions ( |S0D i n or der to avoid 
the central singularity (big bang and big crunch). The term in square brackets depending on x in ( |l02| ) is bounded 
between (as x — ► 0) and ir/2 (as x — ► 2), also, we have c(^ — t ^) — + 2/ (3^/2/7) as k — > 0. On the other hand, the big 
crunch time is ct bc = ct bb + 27rp/n 3 / 2 . Therefore, conditions ( p0[) can only be violated if both k and p diverge but 
complying with < k/ p < 2, a situation that can be avoided if both p and k or pi and ^IZi are finite. As shown by 
figures 0c, Rd and Qe, conditions ( |80| ) are satisfied for reasonable initial value functions, even in the cases displayed 
in fr|d and [fle in which shell crossing occur. In general, a 'good' selection of functions M, K = —E, t bb leads to 'good' 
forms of ^TZi, pi, Yi and viceverza. 

As with parabolic and hyperbolic solutions, initial conditions with density and/or 3-curvature voids lead, in general, 
to shell crossings but a regular evolution with these type of initial conditions is possible in the domain t > U . In order 
to illustrate this point, we consider T given by ( f75| ) and assume that pi > and ^TZi > are everywhere finite (in 
order to comply with (p0|)). Bearing in mind that: 3P(Qi — Q) > — 2 and 1 — P/Pi > for r\>r\i (or t > t{), a 
sufficient condition for T > is given by 



A 



: A 



(to) 



1 

< - 
~ 2 



3 A (*) 
2 1 



(104) 



This condition is illustrated graphicaly by figure ^c, leading to the model displayed in figure |TT[b, an elliptic LTB 
model with an initial density void that has a regular evolution for t > U. 
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D. Simmultaneous big bang. 



An important special class of initial conditions is that characterized by a simmultaneous big bang: ie y = taking 
place in a singular hypersurface marked by t = constant. The conditions for a symmultaneous big bang readily follow 
by setting t' bb = in (|85|), (93) and (103). Clearly, a symmultaneous big bang is incompatible with parabolic solutions, 

since t' = in (p5h implies A < f n ' > = 0, but this case is the FLRW limit (see section VE). From (|93|) and (103), a 



symmultaneous big bang for hyperbolic and elliptic solutions implies the vanishing of the term in square brackets in 
the right hand sides of (|8(]) and (^), that is 



(m) 



(105) 



so that r becomes 



1 +3 Ma, 



(m) 



B A 



(106) 



The conditions for avoiding shell crossings depend on the forms of A and B given by ( p6| ) and (gj) . We assume that 
Pi > and ^TZi > (elliptic) and ( 3 '1Zi < (hyperbolic) are all finite and examine the hyperbolic and elliptic cases 
below: 



• For hyperbolic solutions, both A and B are monotonously decreasing with < A < 1/3 and — 1/2 < B < 0, 
while A — > 1/3, B — > as 77 — > and A — > 0, £> — > — 1/2 as 77 — > 00. For the asymptotic limits, ijkO and r\ — > 00, 

we have respectively T w 1 + A,- m ' ) a nd r w 1 + (3/2) A^ . Since A4 and i?i have also opposite signs, these 
asymptotic limits, together with (105), lead to the following necessary conditions for avoiding shell crossings 



A (fe) > 



2 
3' 



AM A (fe) < 0, 



(107) 



The functions /i, k, A? , A^ are related by ( |92| ) and (105), which can be expressed as 



c 3/2 _ ^2 + x ^fx V2 + x — arccosh (1 + x)\ ^(m) 



c 3 / 2 - |a/2 + x [ yfx^2 + x- arccosh (1 + x)] 



(108) 



2 

1 + - 

x 



1/2 



arccosh (1 + x) 



(109) 



where < £j is the constant value of i bb and a; = |k|//x. 

• For elliptic solutions both A, B are non-negative, with ^4^1/3 and B — ► as 77 — > 0, whil e _B — > (3/2)A — * cxj 
as 77 — » 27r. By looking at the asymptotic behavior as 77 — * and 77 — > 27r and considering ( 105 ), we obtain the 
necessary conditions 



A M > _^ A (m) _ 3 A ( fc) ^ ^ A (m) A (fe) > Q) 



(110) 



As with the hyperbolic case, for t bb = t$ < U a constant, and x — n/fi, the functions /1, k, A^ m \ A| fc ) follow 
from (|0|) and (flOo) 



|(1-£C) - V2" 



3/2 _ |^/2 — x [ arccos (1 — x) — \Jx \J2 — x\ 



(111) 



arccos ( 1 



1/2" 



(112) 
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Since both sets of conditions (107) and (|110|) can be o btained by ap plyi ng (|L05[) to the no-shell-crossing conditions of 
the general cases with t' ^ (ie (po|), (P9|) and (|l00| )), conditions ( 107 ) and ( |lld| ) are necessary and sufficient. 

A simmultaneous big bang greatly restricts the forms of p and k. Notice that it is not possible to obtain a closed 
form k = k{ia) from ( 109 ) or ( |l!2| ). Given a choice of pi leading to p and A^ from (p3|), ( |35| ) and (^l|), the constraints 
(109) or ( |ll2| ) must be solved numericaly in order to determine a very specific and unique form of k for a given 
Having obtained k numericaly, we can use it to construct x = n/p so that A^ can be computed from (108) or (111). 
The remaining functions, such as K and ^TZt follow from ([l3|). This numerical evaluation was done in ploting k and 
A[ k) in fi gures E2I. A simmultaneous big bang leads to a much simpler expression for T, and correspondingly, to a 



significantly more lenient no-shell-crossing conditions (|107|) and (110). Since (107) and (110) are not incompatible with 



positive A or A - fe ' , these conditions can accomodate, besides lumps, some combinations of density and 3-curvature 
voids. Unfortunately, this particular class of initial conditions has a complicated description with the initial value 
functions that we have used. However, a simmultaneous big bang is a restrictive type of initial conditions, like for 
example to assume parabolic evolution (^IZi = 0), or zero density or 3-curvature contrasts (A^ = or A^ = 0). 
We tend to conclude that the small possibility of a dust configuration exactly attaining these very special sets of intial 
conditions makes them somehow unrealistic. Graphical examples of initial value functions and the function y(t. r) for 
LTB models with a simmultaneous big bang and hyperbolic and elliptic dynamics are displayed in figures |l2j. The 
no-shell-crossing conditions obtained and discussed in this section are summarized in Table 1. 



VIII. LTB MODELS 'A LA CARTE' 



The functions Y t , pi and ^IZi play the same role as the traditional variables M, E, t bb , that is, basic "building 
blocks" whose specification fuly determines any given LTB model. It is then entirely equivalent to work with the 
functions M, E, t bb and the solutions of (^) (as is usualy done) or to use the initial value functions Yi, pi, ^IZi 
together with the solutions of (p0|). Given a choice of old variables M, E — —K, t bb , we can find p, Wlli, Yi from 
equations (fuf), dl7j), (|S4|), ( ]92] ) and (102). For either set of variables it is always possible to use one of the functions in 
order to fix the radial coordinate, leaving the other two functions as free parameters. If the traditional variables are 
used in this process, then p and pi can be computed from (^) and (|l6| ) once Y and Y' are known from the integration 
of (JH) . The alternative approach we suggest is to fix the radial coordinate with Yi as in section V C , so that pi and 

^TZi become the free parameters to be used for computing u, k and the contrast functions, A^ , A-^ by means of 
equation s (p5| ) to (fl5|). Next we can obtain T defined in ( p2|) from the solutions of ( p0|) (see section VII ), leading to 
p (from ([2lD) and the remaining geometric and kinematic quantities provided in the Appendix: ^TZ, a, 0, etc. The 
old functions M and E = —K can be determined from (|35|) and ( J3q), w hile t bb can be obtained by setting y = in 
(see equations (|84|), ( 



equations (25) to ( 



and (102) in section VII). 



In the remaining of this section we provide a simple and practical, step-by-step, guideline on how to construct a 
given LTB model in terms of mathematicaly simple ansatzes for the new variables.^. We also show how the character- 
iztion of initial conditions as lumps or voids, as well as regularity conditions and relevant quantities (like y, Y, F, p, 
etc ), can all be illustrated graphicaly by plots that can be easily generated by a standard Computer Algebra System 
(like Maple or Mathematica) . 



Step (1). Prescribe the initial value functions. Consider the following simple ansatzes 

_ a 2 + cup 

Pi - Pi 1 + p ' 



(113) 



(3)ft. =(3)^^ 



b 2 + hf 2 



1 + P 



(114) 



Yi — Si f, 



(115) 



2 These ansatzes are (obviously) not the most general initial conditions and are not based on any realistic observational data, 
we suggest them simply for illustrating the utility of the new variables. The curves plotted in all figures (except and ^) were 
obtained from expressions derived in this section and generated with Maple. 
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where 

• Si is a constant characteristic length and pi, ( s 'TZi are suitable constant reference values of pi, ^TZi 



• / = f(r) is a C° non-negative and adimensional function (see section VC) that satisfies: 

— f(r c ) — 0, where r — r c marks a SC, f(r) > for all r ^= r c . 

— f — > oo as r — > r max , with r max defined by (|49|), 

— If there exist a "turning value" r = r* defined by fir*) — 0, all radial derivatives of pi, ^TZi, p, k, K, M 
vanish at r = r* . 

and the adimensional real constants a\, a 2 , &i, b 2 determine the following generic properties. 

• Central and asymptotic values 

The central values Pi(r c ), ^TZi(r c ) can be related to a 2 , 62 by 

Pi{r c ) = a 2 p l , (3) ^(r c ) = b 2 ^%, (116) 
while in the case when there is a maximal value r max , the asymptotic values /?i(r max ), ' 3 ^7?.i(r max ) relate to 01, 61 

by 

Pi(r max ) = flip, (3) 7^(r max ) = fci^, (117) 

Hence, we can consider the constants 01, 02, 61, 62 as comparative adimensional ratios of special values of pi and 
^TZi with respect to the reference values pi and ^7Zi, for example, values associated to a FLRW background 
(see further below). Numerical values of these constants used in all figures follow this interpretation, hence the 
curves display comparative values of the plotted quantities, say pi, ^TZi, etc, with respect to specific reference 
values (for example, with respect to central values if a\ = b\ = 1). 

• Type of dynamics 

— Parabolic solutions: h = b 2 = 0, so that (3) ft 4 = ( (3) ft 4 ) = K = K = 0. 

— If bi, 62 are both positive, the sign of ( 3 '1Zi is given by the sign of ^TZi. We have elliptic or hyperbolic 
dynamics if (3) 7^ > or (3) 7^ < 0. 

— Mixed dynamics: If b\ and b 2 have different signs, ( (3) ^) and K might pass from positive to negative or 
viceverza (see figure |^b) . 

• Lumps and voids 

Since a 2 / a\ and 62/61 are the ratios of central to asymptotic values of pi and ^TZi, the characterization of ( |113] ) 
and (114) as lumps or voids follows from these constants. 

< a\ < a 2 , density lumps, Q < < a%, density voids, 

< b\ < b 2l curvature lumps, < b 2 < b%, curvature voids, (H8) 

If Pi(r c ) > and ^72-j(r c ) ^ 0, we can identify pi — Pi{r c ) and ^TZi 7Zi(r c ) by setting a 2 = b 2 = 1. In 
this case, we have lumps if a! < 1, b\ < 1 and voids if <zi > 1, b\ > 1. Of course, pi and ^'TZj can also be 
identified with asymptotic values or with values associated with a FLRW background. For example, in a present 
day cosmological context (ti = to, present cosmic time), we can identify pi = po and ^TZi — ^7Zo as constant 
values associated with a FLRW background at t = to. If we set a% = b\ = 1, we can identify these constant 
parameters in terms of the observational parameters fio, Hq of the FLRW background by 



po = B$ti , l^TZo = i? 2 (n -l) 



8ttG 

"3^ u ~ " u "" u ' 6 



The values of po and depend on r (from (113) and (114)), but the other parameters a 2 , b 2 provide the 

ratio of the central and background values: 
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^Wc) = ^-Poa 2 = a 2 Httt , ~ (3) ft (r c ) = ~( 3 )ft &2 = b 2 H 2 (n -l) (119) 

From this example, we clearly have a density lump/void if a 2 is greater/smaller than 1 and a 3-curvature 
lump/void if b 2 is greater/smaller than 1. In all figures we can assume that the numerical values of oi, a 2 and 
b\, b 2 can be interpreted as comparative adimensional ratios with respect to these FLRW reference values. If 
desired, units can be selected so that pi and ^TZi, as well as terms like AirG/c 4 can be normalized to unity. 



Step (2). Compute auxiliary variables: p, K, the volume averages and contrast functions. We need to evaluate the 



integrals (|35| ) and (36) for free functions pi,^ TZt, Yi given by (113), (114) and (115). Let M and K, be defined as 

.3 f (a2 + axf 2 )f 2 fdr 



M = 3 ^ St 



1 + /2 



a i f +3 (a 2 — ai) (/ — arctan /) 



(120) 



/) 



From @, (|§ and we obtain 



3c 4 1 



4nGpi M 



4ttG 



(pi 



(121) 



(122) 



< 3 >7^ AC 
6 / 3 



7< (3) ^>, 



(123) 



A 



(m) 



(og ~ Qi) [/ 3 + (1 + P) (/ - arctan /)] 
(1 + P)M 



(124) 



( fc) = (h - h) [f + (1 + P) (f - arctan /)] 

(i + / 2 )/c 



(125) 



Since /i is simply the average (pi) in units cm -2 , while k and (^TZi) differ by a factor of 1/6, all plots displaying (pi) 
and (^TZi) are qualitatively equivalent to plots displaying /i and k. Numerical factors like or 4TrGpi/c 4 can 

always be absorbed into the numerical values of the constants oi, a 2l b±, b 2 . 

Step (3). Evaluate bang times. 



For parabolic dynamics, the bang time t bb can be computed directly from inserting p given by ( 122 ) into (|84|). 
For hyperbolic and elliptic dynamics with t' bb ^ the bang time, t bb) is furnished by ( |92| ) an d (102 ) in term s of 
k and x. The variable x = n/ p (or x — \k\/p in hyperbolic dynamics) follows directly from (122) and (123) as 



c 4 W% K 



8nG pi 



M' 



(126) 



wher e we must co nside r [^T^J for the hyp erbolic case. The bang time t bb follows directly from inserting k from 
(123) and x from (126) into ( |92| ) and ( |l02| ). This has been done in the plots displayed in figures 0. 



• For the particular case with a si mm ultane ous big bang, the constraints (109) and ( |112| ) cannot be satisfied by 
p and k both having the forms (122) and ( |123| ). However, we can choose either one of these functions and find 
the other one from (10S) to (112). We can choose pt, p, A^" 1 "* from (113), ( 122 ) and (124) and then find k — px 



and A^ fc ' by numericaly solving two of the equations (|108|) to (|112|) . This is the way the plots of figure |12j were 
done. 
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Step (4). Verify fulfilment of regularity conditions 



The no-shell-crossing con di tion s (|S2|), ( |90| ) and (p9|), involving only /i, A\ m ' and Af , can be tested and illus- 
trated directly with (122)- (125) (see Table 1 and figures 0, and 0). However, there is n o simple prescription 
on how to choose oi, <Z2) 6i, 62 in order to comply with the no-shell-crossing condition ( |100| ) of elliptic solutions. 
In this case, it is convenient to expr ess t he functions Pi, Qi in terms of x = k/ jj, by means of Jtt|) , (|78| ) and then 
in terms of ai, a2, &i, bi and / by (12E). As shown by figures [9]a, 0b and [lO|c-e, it is not difficult to select the 
appropriate values of these parameters in order to comply with this regularity condition. 

Regularity conditions, such as (l48|) and (El), take the form 



- = hf 2 + 3(6a-&i) 1 



arctan / 

7 



< 



( 3 )7L S? ' 



(127) 



0. Condition (127) places strong restrictions on 



where the equality holds for every turning value such that /' = 
the values of ^IZi and ^TZi in models with elliptic dynamics. 

Step (5). Choose a specific model. A wide variety of models with specific initial conditions can be constructed with 
(113), (114) and (|115|) by selecting suitable parameters pi, ^Hi, Si and constants a±, 61, b 2 a ssoc iat ed w ith 



comp arative ratios with respect to r eference values of pi and ^TZi (central or FLRW values, as in (116), ( 117 ) or 
(119)), or to lumps/voids (from (118)), or to a given type of dynamics, while, following the discussion of section VI-C, 
simple forms of / can be selected for zero, one or two SC's. Ommiting the cases with mixed dynamics we have the 
following possibilities of fuly regular LTB models: 



Central and maximal values of r are r c — 
or 0b. Bang and shell crossing times are as 



Parabolic solutions with one SC: 61 = 62 = ^TZi = 0, / = tan r. 
and r max = tt/2. Initial density profiles are as shown in figures 0a 
shown in figures 0a and 0b. See figures |a, |b, ^a, |To|b and |ll|a. 

Hyperbolic solutions with one SC: ^7Zi < 0, / = tan r. Central and maximal values of r are as in the parabolic 
case above. Initial density profiles are as shown in figures 0a or 0b. Initial 3-curvature profiles are as shown in 
figures 0a or 0d. Bang and shell crossing times are as shown in figures 0a and 0b. See figures 0a, 0d, [T^a, [iT)| b 
and [Tl|a. 

Elliptic solutions with one SC's: ^TZi > 0, / = tan r. Central and maximal values of r are the same as the 
parabolic and hyperbolic cases a bove . Initial density profiles are as shown in figures 0a or 0b. There are no 
turning values r* , but (f7l|), and (127) must hold. Necessary and sufficient conditions for this are 



0, 



so that: 



K 



1 as 



vr/2 



62 5 



(128) 



therefore, the initial 3-curvature profiles must be as shown in figure [|a, clearly indicating that only curvature 
lumps are possible since ( 3 '7Zi = 2/[Sf (1 + tan 2 (r))] > and so, ^TZ'i < for all r. The bang and shell crossing 
times are analogous to those displayed in figures 0a and 0b. A model with these characteristics but with t' = 
was examined in J5j. See figures 0e, 0f, 0b, ^ and fn|. 



Elliptic solutions with two SC's: ^TZi > 0, / = sin r. Central values of r are r Cl — and r C2 = n. There are 
no maximal values r max and there is a turning value r* = tt/2. Initial density profiles mus t be as shown in figures 
0c or 0d. Initial 3-curvature profiles must be as shown in figures 0b or 0c. Condition (127) must be satisfied 
with the equality holding for the turning value r* 
bi, £>2, Si, ^iZi are linked by the constraint 



7r/2 so that /* = sin(7r/2) = 1. Hence the parameters 



#(71-/2) = 1 



h + 3(62-61) (l 



(129) 



See figures 0c, |d, 0a, 0c, |l0|c, [Tcjd and (□}). 



Since condition ( |12 
turning points f'(r 



cannot hold if K < (parabolic and hyperbolic dynamics), surface layers necessarily emerge at 
= for parabolic and hyperbolic solutions with two SC's. However, these configurations can 
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still be constructed with the expressions derived in this section. 



Step (6). Evaluate and plot time dependent quantities. 

O nce al l pre vious step s have been acco mplished, we have all r dependent quantities in fuly deter mined 

by Jll^)-( pl5|) and (|l2^)-(|T2^) and ( |l26| ), so that all relevant regularity conditions, such as <$^), © , ©, ( |l00| ) 
and (127) can be tested (see Table 1). If this test is passed, we have fuly determined initial conditions for a regular 
evolution. We are ready to derive expressions for y (or Y from (|l8|)) for all t beyond the initial hypersurface %. For 
parabolic dynamics, we ha ve y — y(t, r) from (25). For hyperbolic and elliptic dynamics, we have either an implicit 
solution t = t(r,y) (as in (g6^, ( f27| ) and (|3l|)) or parametric solutions [t — (r,T)),y = yjr, 77)] (as in j2g|), (p9|), f^ ) 
and (|33])). The function T can be obtained as T(t, r) in the parabolic case by inserting ( p5| ) into (|74[), then p — p(t, r) 
follows from (pl|). For hyperbolic or elliptic dynamics we have V = T(r, 77) from ([75]) and p — p(r, 77) also follows from 
(2l|). All o ther time-d epen dent quantities, such as ^TZ, a, O and other invariant scalars given by equations (|A1|), 
( A2|), (A3), ( |A5| ) and (A4), can be transformed as functions of (i, r), for parabolic dynamics, and as functions of 
(r, y) or (r, 77), for hyperbolic and elliptic cases, by eliminating gradients like Y' and Y 7 in terms of T from ( p2[ ) and 
then insterting the previously known T(i, r) or T(r, 77). For hyperbolic and elliptic dynamics, we cannot obtain closed 
functional forms for the time dependent quantities in terms of (i, r). However, we can obtain 3d plots of y (and thus 
of all these expressions) in terms of (t,r), since they are all real functions defined over a domain contained in the 
plane (t,r), hence their image can be represented as parametric 3d plots of the form 



V = [r, t(r,r)), y{t,rf)], 



r 



[r,t(r,r)),T(t,r])], p = [r, t(r,rf), p{t,rf)), 



etc. 



where T(r, 77) is given by (|75[), while p(r,rj) follows from ( |2l| ) once y(r,r]) and T(r, 77) are known (we are assuming, of 
course, that all r-dependent functions, such as pi have already been determined by steps 1 to 5). These 3d plots can 
be easily generated with any general use Computer Algebra System, like Maple or Mathematica. The corresponding 
commands for plotting ?/(i, r) are: 



plot3d([r, t(r,eta), y(r,eta)], r = O..Pi, eta = 0..2*Pi); 



Maple, 



ParametricPlot3D[{r, t[r, eta], y[r, eta]}, {r, 0, 7r}, {eta, 0, 2n}] Mathematica, (130) 

where the parametric functions £(r, 77) , y(r, 77) must be previously defined in each package and we considered in the 
example above the elliptic case with / — sin r. The plotting commands for other quantities follow the same syntaxis. 
The reader familiar with either Maple or Mathematica (or other package) can, of course, embellish the plots by adding 
a wide variety of plot options concerning axes, fonts, labels, colors, etc. All curves and surfaces displayed in the figures 
were drawn with Maple, with those surfaces in figures 0, |, and [ll] plotted using commands like that shown in 
( |130| ). For the case of a simmultaneous big bang, we can still obtain simmilar 3d plots, but t(r, 77) and y(r, 77) need to 
be constructed with numerical values of either one of p, or k. 
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TABLE OF NO-SHELL-CROSSING CONDITIONS. 



In all cases pt and ^TZi must be finite. Also: p; > and: ^TZi > (elliptic), ^IZi < (hyperbolic) 



Parabolic 
Hyperbolic 
Elliptic 



-1 < A (m) < 



-1 < A, (m) < 0, 



-1 < A, (m) < 0, 
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TABLE I. Summary of no-shell-crossing conditions discussed in section VII 
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FIGURES. 




FIG. 1. Comoving layers in elliptic solutions. This figure displays y(t,r) as a parametric 3-d plot [r, t(r, 77), y(r, 77)] (see 
equation (130)) for a regular elliptic model with two SC's. We have considered as initial hypersurface t% = and all initial 
value functions were obtained from section VIII with / = sin r and ai = 1, ai = 2, bi = 1, 62 = 3. Notice how the plane y = 1 
coincides with t = £j = for all layers in the expanding phase, but in the collapsing phase the locuus y = 1 does not correspond 
to a single hypersurface of constant t. The figure also illustrates how the layers in the expanding and collapsing phases in 
elliptic dynamics are (in general) not time-symmetric with respect to any hypersurface of constant t. They are symmetric with 
respect to the hypersurface ij = -K where y = ( hyp ersurf aces of constant 7/ do not coincide, in general, with those of constant 
t). These features also occur for other models with elliptic dynamics, such as models with one SC (see figure 2e) and models 
with a simmultaneous big bang (figure 12c). 



27 




FIG. 2. Comoving layers in terms of Y(t, r) and y(t, r). These figures examine the difference between y and Y by looking at 
the functions y(t,r) and Y(t,r) obtained as parametric 3-d plots for three LTB models: one with hyperbolic dynamics and one 
SC (figures (a) and (b)), one with elliptic dynamics and two SC's (figures (c) and (d)) and one with elliptic dynamics and one 
SC (figures (e) and (f)). Figures (a), (c) and (e) depict y(t,r), while Y(t,r) is shown in figures (b), (d) and (f). Both y and 
Y have been plotted as suggested in (130), with initial value functions given in section VIII having the following parameters: 
figures (a) and (b) [oi = 1, 02 = 3, b\ = 1, 62 = 4, / = tan r], figures (c) and (d) [ai = 1, 02 = 2, b\ = 1, hi = 3, / = sin r], 
figures (e) and (f) [01 = 1, 02 = 15, 61 = 0, 62 = 20, / = tan r]. The difference between y and Y emerges by comparing (a) vs 
(b), (c) vs (d) and (e) vs (f). In figures (b) and (f) we have Y = for the central layer marked by r = 0, while figures (a) and 
(e) clearly illustrate the fact that y > for this layer. The same happens in figures (c) and (d), but now also along the second 
SC at r = 7r. In figure (f) the function Y diverges along all hypersurfaces of constant t as r tends to the maximal value r = 7r/2 
(since / = tan r — > 00). For the same model, as shown by figure (e), there are hypersurfaces of constant t, marked by t < 1 , for 
which y remains everywhere bounded as r — > 7r/2. All these features reinforce the interpretation of y as a local "scale factor" 
for a given layer, in contrast with the curvature radius Y. The grids in the figure correspond to hypersurfaces of constant r 
(comoving layers) and of constant ij. Notice how the latter do not coincide (in general) with hypersurfaces of constant t. The 
initial hypersurface, marked by t = U — 0, is shown as a thick line in figures (a), (c) and (e). Notice that y — 1 for all layers at 
t = 0. The thick curves shown in figures (b), (d) and (f) are the "profile" of Yi — Y(U,r) that corresponds to / in each case. 
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FIG. 3. Initial density lumps and voids. Initial value functions pi (thick curve), (pi) (dotted curve), M and A\ m> (marked 
by the symbol A) are depicted for a hypersurface % with one SC (figures (a) and (b)) and two SC's (figures (c) and (d)). 
Figures (a) and (c) depict initial density lumps in which pi is a local maximum at the SC's, while figures (b) and (d) correspond 
to voids in which pi is a local minimum at the SC's. These functions were obtained from section VIII with parameters: 
[ai = 0.1, ci2 = 1.5] for the lumps ((a) and (c)) and [oi = 1.8, ai = 0.3] for the voids ((b) and (d)), with / = tan r in (a) and 
(b) and sin r in (c) and (d). Notice that (pi) > pi for lumps ((a) and (c)) and (pi) < pi for voids ((a) and (c)), hence A' m ' is 
negative for lumps and positive for voids, clearly illustrating how the sign of this contrast function characterizes the nature of 
the inhomogeneity in density along an initial hypersurface %. In all cases M vanishes at the centers and M' = at the turning 
value r — r* — n/2 (in (c) and (d)), showing the same qualitative behavior for lumps or voids. Since the type of dynamics 
(parabolic, elliptic or hyperbolic) depends on ^IZi and not on pi, the initial density profiles in (a) and (b) can correspond 
to any dynamics (see figures 4a, 5a and 5b). Regular LTB models (without surface layers) having the initial density profiles 
shown (c) and (d) must have the 3-curvature profiles with two SC's shown in figures 4b and 4c. 
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FIG. 4. Initial lumps and voids of positive 3-curvature. The figures display initial value functions ^IZi (thick curves), 

{ (3) Ki) (dotted curves), K and A^ (marked by the symbol A) associated with an initial 3-curvature profile with ^IZi > 0, 
hence the dynamics is necessarily elliptic. Lumps or voids of the 3-curvature correspond to local maxima or minima of ^TZt at 
the SC's. Figure (a) shows the case of a 3-curvature lump with one SC (/ = tan r), while (b) and (c) respectively depict lumps 
and voids with two SC's (/ = sin r). Models with elliptic dynamics must comply with regularity conditions (48) and (54). For 
the functions of section VIII these conditions are equivalent to (127), hence the constant parameters bi, hi, Si and ^TZi must 
satisfy condition (128) for the case in (a) and condition (129) for the cases in (b) and (c). Notice that K — > 1 as r — > 7r/2 in 
(a) and K(ir/2) = 1 in (b) and (c). The contrast function is positive/negative for voids/lumps, while K vanishes at the 
SC's and has the same qualitative behavior for lumps and voids. The initial density profiles that match the configuration in 
(a) is either one of those depicted by figures 3a or 3b, while for the cases in (b) and (c) it would be either one shown in figures 
3c or 3d. 
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FIG. 5. Initial lumps and voids of negative 3-curvature. These figures show ^IZi (thick curves), (^IZi) (dotted curves), 

K and A\ k) (marked by the symbol A), defined by (24), (36), (40) and (42), for the case of negative 3-curvature ^IZi < 0, 
corresponding to hyperbolic dynamics. Figure (a) displays a lump and figure (b) a void. Notice that now lumps/voids 
are characterized by local minima/maxima of ^IZi at the SC (the pictures depict only the case with one SC). The plots 
were obtained with equations (114), (123) and (125), with / = tan r, ^TZi = — 1, Si = 1, bi = 0.1, &2 = 1.5 in (a) and 
bi = 1.5, &2 = 0.1 in (b). The corresponding initial density profile for these figures would be either one of those shown in figures 
3a or 3b. 




FIG. 6. Other initial conditions . Figure (a) depicts the case where pi has many local maxima and minima, so that p\ changes 
sign many times (ie "ripples"), then the contrast function changes sign at each turning value. The function used in this 

example is pi — pi(l + sin r/r) with pi = 1 and there is only one SC. From (48), (54) and the no-shell-crossing conditions 
derived in section VII and summarized in Table 1, such configurations must have elliptic dynamics (with K = 1 at each turning 
value) and might be very prone to develop shell crossings. Figure (b) displays ^IZi (thick curve), { {3) TZi) (dotted curve), K 
and A\ h \ defined by (24), (36), (40) and (42), for a 3-curvature lump in which an initialy positive ^IZi becomes negative. The 
functions {^IZi) and K also change sign (change from elliptic to hyperbolic dynamics). Notice that the passage from elliptic 
to hyperbolic dynamics does not coincide with the zero of ^TZi, hence there is a range of r in which ^TZi < but (^IZi) and 
K are still positive. The plots were obtained with equations (114), (123) and (125) with f = r, ^TZi = 1, Si — 1, bi = —0.4 
and &2 = 1.5. Since (^IZi) and K vanish at the value of r marking the transition from elliptic to hyperbolic dynamics, we 
have (because of (43)) A' fc ' — » — oo for this value of r. 
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FIG. 7. Bang and shell crossing times. The figures display plots of the bang time function t bb (r), equivalent to displaying 
the locus of the central singularity (y = 0) in the coordinate patch (t,r). The locus of the shell crossing (T = 0) singularity is 
also displayed in those cases that exhibit this feature (solid curves in figures (b), (d) and (e)). Figures (a) and (c) correspond 
to initial density and 3-curvature lumps with hyperbolic dynamics and one SC (figure (a)) and elliptic dynamics and two 
centers (figure (c)). Configurations with initial voids are shown in figures (b) (with hyperbolic dynamics), (d) and (e) (elliptic 
dynamics, expanding and collapsing phases). All the plots were made with initial value functions obtained from section VIII, 
with / = tan r and / = sin r for hyperbolic and elliptic dynamics. Lumps with "small contrasts" are characterized by 
oi = 6i = 1, a2 = 1.2, 62 = 1.3, while "large contrasts" correspond to ai = 61 = 1, 02 = 200, 62 = 300. The same values 
characterize voids, exchanging 01, bi for 02, 62- Notice how configurations with small contrasts are "older", in the sense of 
larger time lapse between t bb and U = 0. Also, this "age" varies with r due to the dependence of t bb on r, this variation is 
larger for large contrasts and hyperbolic dynamics (figures (a) and (b)), while for small contrasts and elliptic dynamics it is 
almost insignificant. In figures (b), (d) and (e), the curve of t bb is marked by crosses while those of shell crossings are solid 
lines. Comoving observers (r = const.) in all these cases cannot avoid the shell crossings, but in the hyperbolic case (figure 
(b)) this singularity is confined to the past of the initial hypersurface U = 0. As shown by figures (d) and (e), the shell crossing 
singularity might exhibit two branches for elliptic dynamics, unless we choose initial conditions satisfying (104) as shown in 
figures 9c and lib. The locii of singularities in cases with parabolic dynamics is qualitatively analogous to those of hyperbolic 
dynamics shown in (a) and (b). 
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FIG. 8. Ill-suited initial value functions . The figure displays the bang time t bb for "bad" choices of the initial value function 
pi. The initial hypersurface % is marked by t = 0, we have chosen Yi oc r and ^IZi = (parabolic dynamics), so that t bb 
is given by (84). The curve "A" correponds to pi oc 1/r, a form that is singular at the SC (r = 0), hence % is also singular 
at r = (it is intersected by t bb at r — 0). The curve "B" corresponds to pi oc 1/cos r, diverging as r — > 7r/2. In this case 
7i is singular at r = tt/2 (it is intersected by t bb ). For this choice of pi we have a void with infinite contrast and so shell 
crossings necessarily emerge. Notice (from figure 7) that % is perfectly regular for all initial value functions obtained from the 
expressions of section VIII. We need to select truly pathological initial value functions in order to have a singular %. 
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FIG. 9. Testing initial conditions for regularity. A graphical representation is provided in (a) and (b) for the no-shell-crossing 
condition of elliptic dynamics (100). Figure (c) ilustrates the condition (104) allowing for the possibility of initial voids with 
elliptic dynamics evolving free from shell crossings for t>ti- This type of graphical illustration is not necessary for parabolic 
and hyperbolic dynamics, since the no-shell-crossing conditions in these cases (equations (82) and (90), see also Table 1) only 
depend on the signs of pi, ^TZi and the contrast functions and can be verified simply by plotting these functions (as in figures 3, 
4 and 5). However, for elliptic dynamics the verification of regularity conditions like (100) and (104) is not trivial. In figures (a) 
and (b) we have density and 3-curvature lumps characterized by the expressions of section VIII. The corresponding parameters 
are: [ai = b\ = 1, a,2 = 15, 62 = 20, / = sin r] and two SC's in figure (a) and [01 = 1, 61 = 0, 02 = 15, 62 = 20, / = tan r] 
and one SC in figure (b). The curves marked by "A" and "B" respectively represent the terms 27rP4^ m ' 1 — (3/2)A^ fc '] and 
PiQilA 1 ^ - (3/2)A 8 (fc) ] - [A^ m) - Af >] appearing in (100). For this regularity condition to hold we need to verify that, for 
all r, the curve "A" must be larger than that marked by "B" and both curves have to be larger than zero. As shown in 
(a) and (b) this is what happens for the selected initial value functions (these initial conditions are then used in the elliptic 
models depicted in figures 10c, lOd, lOe and lOf). In figure (c), the curves marked by "A", "B" and "C" respectively denote 
the terms 1/2 + (3/2) Af \ A, (m) and Af 5 appearing in (104), calculated for i nitial value functions characterized by / = sin r, 
a\ = 1.3, a% = 1 (density void) and 61 = 1, 62 = 1.4 (3-curvature lump). For (104) to hold, we need (for all r) the curve "C" 
to be larger than "B" and "B" to be larger than "C" , as it happens in the figure. These intial value functions are used in the 
elliptic model depicted in figure lib. In general, one might choose the parameters ai, 02, / characterizing p^, (pi) and A' m ' 
from section VIII, then we can always find, by trial and error, the appropriate parameters 61 and 62 that comply with (100) 
and (104). 
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FIG. 10. LTB models free from shell crossings. The function F (figures (a), (c) and (e)) and the rest mass density p (figures 
(b), (d) and (f)) are depicted as functions of (t,r) by parametric 3-d plots like in (130). The initial value functions have been 
selected so that shell-crossings are absent for all the evolution of the models (see Table 1). Figures (a) and (b) correspond to 
the hyperbolic model characterized by oi = bi = 1, 02 = 15, 62 = 20, / = tan r (initial density and 3-curvature lumps). Notice 
in (b) the initial density profile (thick curve in the plotted surface) associated with the initial hypersurface t = U = 0. Regular 
models with parabolic dynamics have the same qualitative features as hyperbolic models. The intial value functions associated 
with the elliptic models in (c), (d), (e) and (f) are those tested in figures 9a and 9b. The models in (c) and (d) have two SC's 
(hence / = sin r), while those of (e) and (f) have one SC (/ = tan r as those in figures 2e and 2f). Notice in (d) and (f) the 
initial density profile pi as the thick curve marked by restrcting the 3-d plot to the initial hypersurface t = U = 0. 
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(a) 
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FIG. 11. Models with initial voids that are regular for t>U. These figures depict the parametric 3-d plot of T(t, r) for mod- 
els with initial density and 3-curvature voids having hyperbolic (a) and elliptic (b) dynamics. The model in (a) is characterized 
by ai = 61 = 1, C12 = 0.45, 62 = 0.5. Notice that Y is negative for earlier times but T = 1 along the initial hypersurface 
t = ti = and r > for all t > U = 0. As long as initial value functions are finite, this type of evolution is possible for 
all models with parabolic and hyperbolic dynamics and initial value functions characterized as voids. The elliptic model in 
(b) corresponds to the initial value functions tested in figure 9c. As with the hyperbolic model in (a), the locuus of the shell 
crossing singularity, V = 0, is confined to times earlier than the initial hypersurface t = ti = (marked as a thick line in the 
3-d surface), likewise we have T > for all t > U = 0. As shown in figures 7d and 7e, the locuus of F — in elliptic models 
with initial voids will exhibit, in general, two branches (as the locuus of y — 0). Condition (104), leading to a single branched 
form for T = 0, is a strong restriction on the initial conditions of elliptic models. We used in (b) the elliptic model with two 
SC's but a very similar plot emerges by considering the case with one SC. 
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FIG. 12. Models with a simmultaneous big bang. We consider models with elliptic (figures (a), (b) and (c)) and hyperbolic 
(figures (d), (e) and (f)) dynamics. We select /i and from equations (122) and (124) and find k and A^ by solving 

numericaly equations (108), (109) for the hyperbolic case and (111), (112) for the elliptic case. The functions fi and A^ m ' 
are characterized by oi = 1, 02 —4 and / = sin r, tan r, respectively, for the elliptic and hyperbolic case (we can assume 
that AirGpi/c 1 = 1). For the numerical solutions we have taken the constant bang time in (109) and (112) to be t'fj = 0.23 
(hyperbolic) and t$ = 0.4 (elliptic), for an initial hypersurface marked by t = U = 1. The values of tffl are not arbitrary and 



depend on the choice of ti, \i and A[ m \ Figures (a) and (b) display the obtained k and A\ K> plotted next to the selected fi and 
A'*™' (both contrast functions have the same sign, as required by (110)). Figure (c) displays the parametric 3-d plot y(t,r) 
that corresponds to the full set of intial value functions. Notice that a simmultaneous big bang does not imply (in general) 
a simmultaneous "big crunch" . Figures (d) and (e) display the corresponding initial value functions for the hyperbolic model 
(notice that the contrast functions have opposite signs, as required by (107)). The parametric 3-d plot y(t,r) for the hyperbolic 
case is shown in (e). 



(fc) 
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APPENDIX: INVARIANTS AND GENERIC FEATURES OF LTB SOLUTIONS. 



• Kinematic invariants . For (Q) and we have u a = u a -b u b — 0, so that the 4-velocity is a geodesic field, while 
the expansion scalar and shear tensor are 



e 



y 2Y 



(Al) 



<j a b = diag [0, — 2(7, ct, u] 



oM Y -- Y - 

3 \ Y Y' 



a ab a ab = 6a 2 , 
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Y 



3 +CT = F' 



(A2) 



• Curvature scalars. . The Ricci scalar, ' 3 - ) 7?., of the 3-dimensional hypersurfaces of constant t (rest frames of 
comoving observers) can be expressed in terms of p, O and a by the equation 



.(3)ft 



8ttG 



1 



e 2 



c 2 (jsy)' 



(A3) 



6 '" 3c 2 r 9 3F 2 r' ' 

The 4-dimensional Ricci scalar, Ricci tensor squared and Riemann tensor squared (the Kretschmann scalar) are 

n ab n ab = K 2 

4 /3m \ 2 5 



_ 8ttG 



r 3 



3 P 



(A4) 



while the electric part of the conformal Weyl tensor, E 10 & = C a bcd,u a u b , is given in terms of the only nonzero 
invariant Weyl curvature scalar, ip2, by 



E a b = diag [0,-2^2,^2,^2], 



^'2 



4ttG 
"3^" 



3 m 
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F 3 



4ttG 
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Regularity and differentiability. 

Spherical symmetry is characterized by the group SO (3) acting transitively along 2-spheres (orbits) with proper 
area 4-7T Y 2 (t, r). The domain of regularity of LTB solutions can be given then as the set of all s paceti me p oints 
marked by (t, r) for which (a) Y > and (b) all curvature scalars (like those in (A3), (A4) and flA5| )) are 
bounded. 

Minimal differentiability requirements for the basic variables can be summarized as follows fl6[| , [H: Y and 
all its time derivatives are continuous, while quantities involving radial derivatives like Y', Y' , Y", etc can be 
piecewise continuous, thus allowing jump discontinuities for a countable number of values r = r\ so that the 
right and left limits as r — > r\ are different but this difference is always finite. These conditions imply (via 
(§)) that the basic variables Y, M, E, t bb must be at least continuous (ie G°), while p and other quantities 
(0, (7, ^TZ, containing terms like Y' or Y' could be piecewise continuous. 

From (0), the weak energy condition leads to 



p > <=> sign(F') = sign(M'), 

while regularity of the metric component g rr requires 

1 + E > 0, 



(A6) 



(A7) 



a condition that is trivialy satisfied for parabolic and hyperbolic (E > 0), but not for elliptic solutions or regions 
(E < 0). The equality in (A7) implies that Y' and E must have a common zero of the same order (see below). 
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Zeroes of Y and Y' 

The functions Y and Y' can vanish, either for a value r — const., or for a coordinate locuus (t, r) with r =/= 
const. The two cases lead to very different situations: 

— Zeroes of Y for r = r c (Symmetry centers, SC). 

These regular worldlincs marked by r — r c are the worldlines of fixed points of SO (3). LTB solutions (like 
all sphericaly symmetric spacetimes) admit zero, one or two SC's. They are characterized by Y(t, r c ) = 
Y(t,r c ) = 0. All of M(r c ), M'(r c ), E(r c ), E'(r c ),t' bb (r c ) vanish at a SC, but the order of the zeroes of 
these functions (power of leading terms) follows from the regularity of each term in (||). Near r = r c we 
have g|, @, @ 

Y*Y%r B ){r-r c ), M « M'"(r c )(r - r c ) 3 , E m E"(r c ){r — r c ) 2 , t bb » t£(r c )(r - r c ) 2 , (A8) 

— Zeroes of Y for non-comoving locii (Central singularities). 



All curvature scalars, such as (|A3j), (|A4|) and flA5| ), diverge as the proper area of the orbits of SO(3) 
vanish, hence the name "central" (associated with singular centers). By analogy with FLRW spacetimes 
and depending on whether dust layers expand or collapse, these singularities are called a "big bang" or a 
"big crunch". These singularities cannot be avoided for any choice of free parameters, thus they are an 
essential feature of LTB models. 

Zero of Y' for r — r* (Turning values of Y 1 and surface layers). 

The extrinsic curvature forms of the 3-dimcnsional submanifolds r = const, must be continuous (at least 
C ) as r increases. If if Y'(t, r*) = for a comoving r* ^ r c , these extrinsic curvature forms are continuous 
as r — > r* as long as the following regularity conditions hold Q] , |l6| , , |Q 

Y'(t,r*) = O M'(r*) = 0, 1 + E{r*) = 0, (A9) 

implying that the common zero, r = r* , of Y', M', 1 + E is of the same order in r — r* . If these conditions 
fail to hold, the extrinsic curvature of hypersurfaces r — const, is discontinuous as r — > r* and we have a 
surface layer at r = r* (see @, @, @, @, @). 

Zeroes of Y 1 for non-comoving locii (Shell crossing singularities). 



From (|J), energy density p diverges (and thus all curvature scalars like (|A3|), ( A4) and (A5) diverge). Since 
g rr vanishes for (in general) Y > 0, proper distances between comoving layers vanish, hence the name "shell 
crossing" . Unlike central singularities, shell crossings are unphysical and are not esential. Conditions for 
having LTB models free from these singularities can be given in terms of M , E, t bb and their gradients (see 
section VII, Table 1 and references ||, @, 0, @). 

"Open" or "closed" models. Hypersurfaces t — const, T, are everywhere orthogonal to the 4- velocity vector 
field u a , hence they can be characterized invariantly as the rest frames of comoving observers associated with 
the metric 

hab = u a u b + g ab (A10) 

All T are Cauchy hypersurfaces in LTB solutions. They are foliated by 2-spheres marked by constant values 
of r. However, the radial coordinate can be arbitrarily rescaled, and so an invariant measure of distances along 
its direction is given by the proper length integral £ = J ^/g rr dr, evaluated for spacelike curves with fixed (9, 4>). 
For a regular and maximaly extended hypersurface, the domain of r is limited, either by SC's (since Y > 0) 
and/or by maximal values r = r max defined by 

f Y 'dr /A , 

1=1 — > oo as r — > r max , (All) 

J VTTe 

The terms "open" and "closed" respectively corresponds to an infinite or finite proper volume of these hyper- 
surfaces, evaluated from (]l|) as 

r Y 2 Y'dr 

v - 4 */m < A12 » 
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where the integral must be evaluated along all the domain of r. The class of homeomorphism of hypersurfaces t = 
const, is related the number of SC's (see section VI-C for more details). From familiarity with FLRW spacetimes, 
one would tend to expect a "natural" association between closed topologies with elliptic (recollapsing) dynamics 
and open topologies with parabolic or hyperbolic dynamics. However, this association need not be mandatory 
(see [|l2| and sections VI and VIII): open topologies might have elliptic dynamics and closed topologies 
might present parabolic or hyperbolic evolution, though in this last case a surface layer appears at r = r* such 
that Y'{t,r*) = (sec |§, §, @ and @). 
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